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1 164 EVELOPMENT 

‘WORDS: Currents (water); Dispersion; Lagrangian functions; Lake 

= Superior; Numerical analysis; Plumes; Remote sensing; Runoff; 

| Sedimentation; Suspended so solids; ; Turbidity; Water pollution; Water pollution | ; is 


ABSTRACT: Currents as a function of wind were calculated by a depth | integrated — 
model of Lake Superior. Using the calculated currents, the development of large ‘ 
i | oo plume was numerically simulated. Good agreement between the simulated 
ie. remote sensing, and ground truth data provided verification of | calculated — 
_ current patterns and much insight into the ultimate fate of pollutants in western Leke 
‘Sener as a function of wind. The long range transport of a conservative pollutant - 
from the Duluth-Superior harbor was simulated for westerly winds and for variable 
winds. For westerly winds, the model predicted that the contaminant would be ~ 
_ Wansported along the Wisconsin shore away from the harbor area. For variable winds, 4 


pollutant remained trapped in the southwest corner of the lake. “ot 


_REFERENCE: J. (Assoc. Scientist, Dept. of Physics, Univ. 

Minnesota at Duluth, Duluth, Minn. 55812), and Sydor, Michael, “Plume * 
- Development Using a Lagrangian Marker Method,” Journal of the Waterway, Port, — 
_ Coastal and Ocean Division, ASCE, Vol. 107, No. WW3, Proc. Paper 16462, August, ym 


1981, pp. 131-148 atensive ate gene 


- 16460 DISTRIBUTION OF CREST-TO-TROUGH WAVE HEIGHTS 
KEY WORDS: Field data; Probability density functions; Probability 
| distribution functions; Rayleigh waves; Troughs; Wave height; Ww ave 
measurement; Wave spectrum © PALE 
Late lor winds <3). be more fully the. 
| The theoretical distribution of wave heights as given by the Rayleigh al “4 


t 


is based on the wave envelope. However, empirical distributions derived from field | 
data utilize a crest-to-trough representation. The latter definition corresponds to an 
ri averaging operation on the envelope which results in an excess of waves with heights 2 | 
near the midrange and in deficiencies at the two extremes in a manner increasing with 7 ; 
| the spectrum band width. The ultimate effect of a deficiency towards the high wave oi 
1 tail is an underestimate in the exceeding probability of larger wave heights. Hence, this _ 
! kind of discrepancy (systematically observed between the Rayleigh law and various ad 
| expleined on this basi, of advection from 
| REFERENCE: Tayfun, | M. Aziz (Assoc. Prof., Civ. Engrg. Dept., Kuwait Univ., P. oO 
7 Box 5969, Kuwait), “Distribution of Crest-to-Trough Wave Heights,” Journal of the 
Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 1 107, No. ww, , Proc. ale 
| 16460, August, 1981, pp. 149-158 
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16461 EXTENDED VELOCITY POTENTIAL WAVE KINEMATICS _ 


KEY WORDS: Boundary conditions; Coastal engineering; Reet — 
Fourier series; Ocean engineering; Stokes law; Velocity; ‘Water depth; Wave 


ht; W Wa at WwW locit 


ABSTRACT: - double Fourier series expansion a of the velocity potential for water 
waves was used to treat nonlinear single harmonic waves and irregular waves measured — 
in the ocean test structure experiment. The (extended) velocity potential | (EXVP) 
satisfies the governing hydrodynamic equations in the body of the water and matches — 
the boundary conditions at the waves surface in the least-squares sense. The waves that 
are treatable with this velocity procedure can have a range of crest values for a 
specified wave height, period and water depth. In contrast, the common Stokes — 
velocity procedures can treat only waves of fixed crest, for a given wave height, period 
and water depth. The water velocities measured for irregular waves in the ocean test 
structure experiment are in good agreement with the velocities predicted by this 
procedure; these correlations are significantly better than velocity correlations with the — >. 


+ 
REFERENCE: Lambrakos, Kostas F. (Sr. pene Specialist, Exxon Production 
Research Co., Houston, Tex. 77001), “Extended Velocity Potential Wave Kinematics,’ 

Journal of the Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 107, ‘No. 


WW3, Proc. Paper 16461. August, 1981, pp. 159- 
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16478 OSCILLATORY ROUGH TURBULENT BOUNDARY I LAY 


WORDS: Boundary “layer; Fluid mechanics; Hydrodynamics; 


Laboratory tests; Mathematical analysis; Mathematical models; Turbulence; - 
boundary layers; Turbulent flow; Water flow + bal 
ABSTRACT: Kajiura’s s layer ‘model (1968) « of an oscillatory rough turbulent layer ‘ 
is simplified to a two-layer model. This model consists of an overlap layer in which the 
_ turbulent viscosity varies linearly with height, and an outer layer in which ne viscosity 
is constant. The analytic expressions are compared with Jonsson’s Test No. 1 (1963, — 
1976, 1980). Reasonable agreement between theory and and experiment is found, thus — 
REFERFNCE: Brevik, (Lect., Division of Port and Ocean’ Baginecring, 
University of Trondheim, Trondheim, Norway), “Oscillatory Rough Turbulent 
Boundary Layers,” Journal of the Waterway, Port, Coastal and Ocean Division, ASCE, 
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USING LAGRANGIAN 4 


‘general current patterns from locally turbulences. 
modeling a particular event, such measurements are costly and difficult to obtain. ae 
Sequential remote sensing data on the transport of a large turbidity plume provides: - 
an alternative to extensive current measurements. In the past, remote sensing 
data has been used to qualitatively verify calculated current patterns in western 

_ Lake Superior for various winds (3). _ In order to ‘more fully utilize the pater, 1. 


of remote sensing data as a means of verification of calculated current patterns, — 


we numerically simulated the development and transport of the ii aa 
plume from the Nemadji River, Wisconsin (Fig. 1), 


ole an accurate numerical solution of the equations which govern the 
1 transport and diffusion of turbidity and a detailed knowledge of all major sources: > 
% and sinks of turbidity. Most numerical solutions of advection suffer from distortion 4 


_ The use of plume simulation as a means of verifying calculated er a 4 


near the edge of the plume where concentration gradients are necessarily large. ~ 
_ In a first attempt at plume simulation, an explicit integration of the advection- 
_ diffusion equation was used. The problem of distortion was alleviated by use 
of upstream concentration gradients | at all points, but the numerical scheme sl 
still suffered from excessive dispersion at the edge e of the plume. 
a to avoid this problem a a Lagrangian marker technique was tried (14). This solution © 
- technique was ideally suited for modeling of a point source of pollution and 
proved of handling concentration gradients without significant 
distortion. 


uel and sinks of turbidity. The major sources es of red clay tatty in n western ; 
Lake Superior are shore erosion, resuspension, and river runoff (16,21). During _ 
ar the 1976 spring runoff period, March 31—April 6, river runoff was the most — 


a Assoc. Scientist, Dept. of Physics, Univ. of Minnesota at Duluth, Duluth, Minn. 55812. 
ee Prof., Dept. of Physics, Univ. of Minnesota at Duluth, Duluth, Minn. 55812. MOS 
oe Note.—Discussion open until January 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, © 
_ ASCE. Manuscript was submitted for review for possible publication on April 15, 1980. 
This paper is part of the Journal of the Waterway, Port, Coastal and Ocean Division, 

Proceedings of the American Society of Civil Engineers, -©ASCE, Vol. 107, No. WW3, » 

August, 1981. ISSN 0148-9895 /81/0003-0131/$01.00. 
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Numerical lake models are nor mally tested Dy Comparing their results with 

— 

— 

| | 


substantial shore erosion, and strong northeasterly winds needed 
generate significant resuspension occurred only on April 3. The honey River Bee ae 


western arm of sah teentiler: As part of our effort for PLUARG (Pollution a 
from Land Use Activities Reference) we monitored the 1976 spring runoff peak — 
of the Nemadji (18). These observations together with a rough approximation ae % 
; ; at es resuspension on April 3 made possible a fairly accurate simulation of the bt 
; ‘major sources of lake turbidity. | The only sink of lake turbidity is settling. 
-* Consecutive Landsat images on | April 5 and April 6 indicated the settling rate oF] 
was about 6% per day during the Tunoff period. This settling rate is 8 consistent — 


a 


Douglas unity Wisconsin od rig) 
ix 
large lakes the movement contaminants is dominated 
4 f advection (4,6). In anumber of model studies, including the present one, horizontal - 
_ diffusion lias had little effect upon model predictions (4,7,11). Since our seme 
_ simulation model incorporates an-accurate, unconditionally stable numerical 
ee ‘solution of the advection diffusion equation with a detailed knowledge of the 
a major sources and sinks of turbidity, its accuracy should therefore depend he 
primarily upon the accuracy of the current patterns it employs. An accurate “ 
simulation of plume development would confirm calculated current patterns. > a 
_ - Numerical simulation of turbidity transport is useful in other ways. Suspended 
i ‘ted clay is one of the major pollutants in western Lake Superior. In addition, “ “ee 
the highly turbid water from the River i is easily detectable in 


PLUME DEVELOPMENT 133 
harbor. Together ‘with the St. Louis River, | the Duluth- -Superior constitutes 
the major United States source of man-made pollution into Lake Superior. — 
the major simulation of plume development yields a better understanding _ 
the processes involved in pollutant transport and provides much ins much insight | into” 
the ultimate fate of pollutants as a function of wind ~~ rs anes ee 
_ The 1976 spring runoff of the Nemadji River was monitored by several groups 
under project PLUARG. High altitude manned overflights by NASA’s Lewis 
_ Research Center, Cleveland, Ohio, provided multispectral remote sensing data 
of western Lake Superior on April 1, 2, 4, 5, 6. The western arm of the lake 
was also visible in consecutive Landsat satellite images on April 5 and 6. = 
i were taken by ship cruises on March 31, and April 1, 4, 5, and 6. a 


Sampling stations were to monitor of clean 
a lake water, St. Louis River efflzent, Nemadji River effluent, and lake water 
alon f; 
g the Wisconsin sh¢ shore affected by resuspension, erosion, and output from ; 
minor rivers. At each station vertical ‘profiles of turbidity, temperature, and 
- conductivity were measured. Samples were taken at the surface, the mid-secchi 
layer, and the bottom ami later analyzed for suspended solids concentration | 
a The data from the sampling cruises are tabulated by Sydor and Oman (18). _ 
As Storz et al. (16) have shown, suspended red clay concentrations in western _ 
Lake Superior and Landsat remote sensing data are extremely well correlated. 
Sydor and Oman (18) were also able to identify various types of particulates, : 
using Landsat data. In correlating their remote sensing data from April 1, 
4, 5, and 6 to measured concentrations of suspended solids, NASA’s 4 
Center, obtained correlation coefficients varying between 0.93-0.99 


| 


station was located near South Superior, Wisconsin, at the upper end of the | 
estuary zone. Measurements of suspended solids, stage, and flow rate were 
taken at least once every day during the entire runoff. The suspended load 
output of the river for the runoff was derived (Fig. 2). The total load output 


for the runoff was estimated at 60,000 metric tons. This figure pertains cer 


to the red clay particles less than 60 y» in diameter. Generally, particles smaller 
than 4 4% comprise the major part of the immediate contaminant output into — 
Currents as a function of winds ° were calculated by a a depth integrated model 
of Lake Superior with quadratic bottom friction following a numerical scheme 
by Leendertse (5). The basic methods and calibration of the model are considered 
: by Diehl et al. GB). In this model, Lake Superior v was approximated wang a 


Fig. At boundaries between the lake ond shore, zero transports were maintained. 
The model was started from rest during a period of light winds on March 29° 


at 2200 hrs CST. Using a 5 min time step, the model was stepped forward v 


—Staggered Grid agree in Lagrangian Marker Particle Model “8 


| actual suspended red clay concentrations could be derived from remote sensing 4 - 
“ 
| 
% 
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“Atmospheric Administration (NOAA) wind data seven stations 
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AG. 4 —Current Pattern for Westerly Winds Predicted by Depth Int integrated sa 
Model: hee may be replaced by The bar 
aquatic can thes ba writen 


distance between the grid point and the stations as weighting factors. The surface q 


a weighted average of the measured wind velocities using the squared reciprocal 
Ss was assumed proportional to the wind speed squared. = model —_ 


pnt during the runoff period, 
7 
. 
| 
| 
&g 
4 
—_ 
7 
it 
| 
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calibrated by varying the wind sentn cocfficient so that the we of predicted i= 
and measured water-level oscillations around the | lake matched. Aircraft overs 


flights, land sightings, and observations of the plumes during a cruises  & 
qualitatively verified the predicted transport directions. 
In the results of the lake model, two distinct current pares 1s appeared in 
Ps extreme western arm of Lake Superior. From March 31-April 2, northwesterly ; 
winds 5 mph-15 mph (2 m/s-7 m/s) produced a counter-clockwise circulation — 
_ pattern in which water came into the extreme western arm of the lake a —s\ ; 7 
the north shore and left along the south shore (Fig. 4). On April 3, a moderately — a a 
strong 25 mph (11 m/s) northeast wind rapidly shifted the current pattern to = 
one in which water came in along both shores and left along the axis of = aq 
lake (Fig. 5). Westerly winds returned on April 4 and by April 5 had again = 
_ produced a 4 counterclockwise circulation pattern. . The general character of these = 
4 circulations was brought out by Diehl et al. (3). 
3 order to model the Nemadji River turbidity plume i in detail, the 3.7 mile - 
x 3.7 mile (6 km x 6 km) grid used in the calculation of transports was subdivided 
into a 0.6 mile x 0.6 mile (1 km x 1 km) space staggered subgrid extending 
from Duluth to Bark Point. A subgrid depth array was estimated from depth 
data read off Coast Guard Chart LS 96. Subgrid current arrays were — 
from bilinear interpolation of the 3.7 mile x 3.7 mile (6 km x 6 km) current 
arrays. The interpolated current arrays were then modified to make them meet 


the continuity equation on a 0.6 mile x 0.6 mile (1 km x | km) basis. os ih 


Since the current model predicted only depth integrated currents, vertical: 
movements and gradients could not be taken into consideration. Since vertical | 
profiles of | concentrations of fine red clay particulates are generally y observed 
to be fairly uniform in the shallow waters where lake turbidity is usually found 
(16), the equations which govern the transport of pollutants could be integrated a 
over the vertical without significant loss of accuracy. The depth | integrated 
concentration of suspended solids C, was then simulated using a two- ‘dimensional - 

‘model. Concentrations: were found later by dividing by the depth. A eS 
_ A Lagrangian marker particle technique was used to simulate turbidity transport — 
(14). In order to derive the basic principles of this method, the three-dimensional 


Sy 
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The diffusivity matrix D D is , assumed to be diagonal. Letting D, denote the 
horizontal diffusivity, Eq. | can be integrated from the bottom z = —H to 
the free surface z using the hydrodynamic the 


‘mathematical definitions of surface and bottom flux: ad 


be 


| 
ar 
> 
q 
¢ 
equation a Conservative substance in an incompressible 
) 
&g 
— 


In this equation, n, S represents" the sum of the diffusive fluxes into and out “— 7 


7 of the water column at the surface and bottom. In a two-dimensional model, — 


Ss can be thought of as a source term. Since vertical diffusion is assumed to a 


_ be rapid enough to maintain uniform concentration profiles, the he following depth — 
_ integrated parameters may be defined: 


edz = =(H me wi) wiles 


sent 


depth averaged concentration; and U and V = the depth- averaged components ' 
of the current. Since h << H, H + h may be replaced H. The horizontal 


-D, 


Using bed ome advection-diffusion equation ca can 


3 an 


6 has the following physical interpretation. The quantity is being added 


at a rate specified by S and is then being moved with velocity V;. This _ 


interpretation forms the basis of the Lagrangian marker particle method. In 
Eulerian methods, the region of interest is divided by a grid and contaminant 


q 


Py 

equation can then be written in 
_ 
_ Following Sklarew (13) the net velocity V,isdefined 
| 

concentration in each Cell is tabulated. At cach time step, material convected =f 

_ from neighboring cells is completely mixed within a given cell. This mixing § 


can lead to false numerical diffusion if an ‘inaccurate er of the ‘finite 4 


_ difference equation is used. Eulerian methods with little numerical diffusion 
. g often suffer from spurious oscillations near the edge of a plume where concentra- 
7 tion gradients are large. The Lagrangian marker particle method solves this 
problem by dividing the contaminant load into a number of packets or marker = 
particles. These particles are introduced into the lake at contaminant sources — a 
at a rate corresponding to contaminant loading and then allowed to move with 
4 In the present model, the region of interest is first divided by a square Eulerian 7 
ns and the depths at the center of each grid square are specified. Locations 
_ of contaminant sources are specified and marker particles are introduced at e 


= these points. Each particle is assigned a size smaller than the Eulerian rid a 


the contribution of a marker particle to a Eulerian grid square is given by 
the product of the area of oveilap and the e depth integrated | concentration | of 


_ The contribution of diffusion to the net transport velocity is then added to. 
currents provided by he hydrodynamic model. Because the grid is staggered - 
in space, whenever C, or H in Eq. 5 is zero, its corresponding derivative is 

bai zero — the quotient is s is set to zero. The marker particles a are then moved 


a and depth- stoma concentrations of the manner particles. In this calculation, 


| 
4 
— 
> 
_ 
= 


_ by the net transport velocity. A counties: ates is used to assure be 
7 particles move with the current. If the position of a a particle is given by Sua: ; ; 
the current V,,, is found by bilinear interpolation of current arrays. An —. 4 


= old 


- This iterative sch scheme has been shown to greatly i improve the accuracy of particle 
on whether they were above « or below background concentrations. Based upon 2 
analogy with an analytic solution of diffusion of a square patch (15), the size 


2 =6D o(t+t,). 


Diag Eq. with respect to time 


of the ‘concentration The change. in size of the particle v was. 
given by a function of C, and the concentration of particle 
4 
In order | to avoid numerical difficulties when A approaches zero in Eq. ll, 
particles were never allowed to get smaller than an eighth of the Eulerian grid 3 
size. If particles became larger than the Eulerian grid size, they were subdivided _ 
_ to make four new particles. This scheme made possible a more accurate 
representation of the dilute fringes of the plume. 
- With the new particle positions and sizes, a new concentration | array is eines = 
and the whole procedure is repeated. ‘This numerical scheme was capable of 
accurately modeling a runoff plume at all stages of development with a limited 
_ number of particles and a limited amount of psig renee time (about 10 min on > 
a CDC CYBER 74 for a one month simulation). 
Observation of individual patches in consecutive Landsat images on April 
. 5 and 6 made possible direct measurements of horizontal diffusivity and settling 
tates (12). Based upon these studies, a a constant horizontal diffusivity of A 
sq ft/sec (1 m’/s) was used in plume simulation. This figure is consistent 
with empirical formulas for horizontal ag) 


i 
particle position Raw is found by Ea. 
} 

> | 

' 
| 


plume of the Nemadji River was numerically simulated. The ‘marker particle 7 
model was started on March 31 at 0000 CST and integrated using a 30 niin 
s step. Marker particles were introduced every time step. The suspended 
: _ load input at the Superior entry followed and measured output of the Nemadji_ — 
for the 1976 spring runoff. Based on studies of the 1975 > spring runoff (21), — 
: 7 a suspended load equal to 5% of the ‘Nemadji output was introduced at the | 
mouth of the Amnicon River. Input prior to March 31 was approximated by 
an initial 3,728 tons at the mouth of the Poplar River, I, 912 tons at the mouth aN 
of the Amnicon, and 4,875 tons at the Superior entry. A settling rate of 6% 
rs per day was allowed during the peak runoff period. Since high turbulence and — a 
resuspension» are associated with strong northeast winds (21), all 


3 
4a 


---- 10 mg/ Simulated 


7.—Simulated Turbidity Line) 10 mg/L 
for April 1, 1976, Compared to Surface Suspended Solids Concentration (Solid Line) ; 
_ Obtained from Manned Overflight Multispectral Data in Correlation with = 


Measurements (Simulated Plume Takes Only Sources at and rs 
into Consideration 


settled inavedel up toa maximum of 50 metric tons per grid square was resuspended = 
_ at noon on April 3. Since a protective ice shelf prevented significant erosion 
= the runoff peak, this scheme of introducing and removing suspended — 
load accounted for most of the major sources and sinks in the western arm 
@ the lake. Some material in the actual plume came from saenapension of 
residual material along the Wisconsin shore. This source was neglected. _ 
.- - predicted plume followed the observed plume reasonably well. From 
March 31 to April 2 the predicted plume (Fig. 7) filled the southern corner 
of the lake, as was anticipated from currents shown in Fig. 5. The NASA 
<n and boat observations from that time showed a well-defined high 
density plume moving south along the shore. On April 4, the model and manned 
overflights clearly showed (Figs. 8 and 9) the ‘Sweeping of the plume straight 
out along the a axis lake the northeast winds on April 3. 
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* PLUME DEVELOPMENT 


_ moving east. The predicted plume for April 6, (Fig. 10) displayed many of 


= observed in the satellite image for that day. The 


od 


» FIG. 8.—Simulated Plume for April 4, 1976 (Broken Line) and Suiléie Suspended 
Solids Line) from M Menned Overtlights, at 5 Concentrations 


a 
--- Simulated sus. sol. uy, 7 


numbers show depth av. sus. sol. 
~ 


FIG. | 9. —High Suspended S Solids ‘Simulated | Plume (Solid | Line) for 20 mg mg/L 
Concentration, Compared with High Density Plume Represented by 1 m Secchi 4 
Transparency Contour (Broken Line) from Manned Overflights in Correlation with = 
Sampling Data (Numbers Represent Averages of ‘Mid- secchi and Bottom 
plume, however, contains a certain amount of material from resuspen- 
sion along the south shore which, as pointed out, could not be simulated for 2 
_ lack of data. The total anspented | load in the Landsat image on April 6, 29, 000 


| 
| 
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tons, agreed well with 27, tons in n the plume. 
~@ number of factors contributed to the differences in the observed and predicted | 
contour lines in Figs. 7-10. First, the initial distribution of turbidity and 
resuspension on April 3 could only be crudely estimated. Second, features in 
“the real plume generated by current patterns smaller than the 3.7-mile (6- km) — 
grid size of the hydrodynamic model could not be accurately predicted. Third, 
a small difference between real and predicted currents would be amplified in 
e time. For instance, a constant 0.03 ft/sec (1 cm/s) difference would shift the 
plume Position 0.5 mile /day (0.9 km/day). 


_ However, despite these Sources: of error, _the turbidity model realistically 


solids as a ‘function of time. This 2 accuracy y verified the general current patterns 
but not the point by point currents the 


- for April 6, 1976, 10 mg/L Concentrati 4 
Compared to 10 mg/L Surface Suspended — Obtained from Landsat Data | for _ 
April 6 in Correlation with Sampling Measure ents 


order to obtain a better understanding of salitines transport in western Lake 
7 Superior, the plume simulation model was also used to predict: (1) The distribution dy 
| accumulated red clay sediment; (2) the distribution of runoff water from 
the St. Louis and Nemadji Rivers; — long-range transport of contaminants ~ 


_ The distribution of sediment accumulation is important because ‘Tesuspension — 


of previously settled material is one of the major causes of poor water quality 
_ in western Lake Superior. Resuspension appears to be largely confined to an 
area along the Wisconsin shore east of the Superior entry. Generally, the 
y near-shore zone up to a depth of 65 ft (20 m) is scoured (1,2) with the exception 
of the area east of Wisconsin Point where thick deposits of red clay are found © 


in shallower waters. The predicted distribution of accumulated sediment on . 
6 is Consistent with these observations. area of maximum 


— 
q 
| 
Suspended Solids April 6,1976 / 
he 
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DEVELOPMENT 


deposition of oe east of Wisconsi Wisconsin Point was the result of an eddy aii 


in the southern corner of the lake. Since the near shore zone is scoured, about 


—— 50% of the runoff sediment would later be resuspended and moved to deeper 
_ In a second application of the model, the discharge from the St. Louis and bs 
Nemadji Rivers was traced in order to simulate the transport of a conservative _ 
pollutant from each of these rivers. The model was modified to keep track 
_ of the fraction of river water in the water column. For | the Nemadji, the input 
- at the south entry followed the river’s measured discharge. The initial distribution 
was approximated by 14,000,000 cu yd (11,000,000 m*) of Nemadji water near 
the south entry. Based mpon the results of a numerical model of the md 


— (19), 60% of the St. Louis River’ s measured discharge was loaded at the north 
entry and the remaining 40% at the south entry. The initial distribution o of St. 
1 Louis | nent consisted sail 41,000,000 cu yd (31,000,000 m yn near | the north 
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entry y and 26,000, 000 | (20,000,000) near | the south e entry. The discharges from 
each of these two rivers during the peak of the 1976 spring runoff were traced 
separately using the seven days of calculated currents. . The results compared — 
reasonably well with Landsat imagery and f field observations. _ Fig. 12 shows 
calculated lines of 10% St. Louis River water and the 5 mg/L suspended solids 
line for simulated variable wind conditions three weeks after the spring runoff 
of Lani! dow tet carty 

In a third application of the model, the long-range transport of contaminants 

- western Lake Superior was simulated by tracing St. Louis and Neen 
River waters for one month. In order to reduce costs, the time € step was increased 


3 

¥ to 2 hr. Based upon past : studies, the settling ‘rate was increased to 15% Qn). 
a cases were run. In one case, the seven days of calculated currents were _ 
recycled, in effect simulating one month of transports which could arise as 

‘the result of the second case , only y current patterns Given, 


| 
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| 
7 «FIG. 11.—Simulated Accumulation of Sediment Distribution from Spring Runoff up 7 


by westerly wind patterns were used. The results for simulated effluent dispersion . 
for both cases are shown in Fig. 13. In the case of variable winds, the runoff 


"remained west of the Brule River. Ite is anticipate’. thet that the” effluent would 
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SIMULATED ST. LOUIS RIVER RUNOFF AFTER 30 DAYS a / 
WINDS 
VAR 


A 
LANDSaT DATA MAY 3, 
TURBIDITY PLUME [3NTU] 
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J FIG. 13.—Simulated Distribution of Runoff from St. Louis One Month after Onset — 
of Runoff: Broken Dotted Line Shows Contour for Variable Winds; Broken Line “ 

Shows Contaminant Distribution for Westerly Winds; and Solid Line SI Shows © Outline : 

of Turbidity Plume for for May 3, 1976, , One Month | after 19; 1976 Runoff 30) 


= 
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‘ 
LAKE SUPERIOR 
x : FIG. 12.—Simulated St. Louis River Effluent and Red Clay Distribution for Variable 
on April 20, 1976, Three Weeks after Nemadji Runoff 
= 
4 
J 
ihe ency contaminants to accumu ate in| western arm ake oupe 


math 
yy hes been inferred from the investigations by Ruschmeyer et al. (10). The apparent - 


_ trapping of runoff water is also evident from Landsat images. In the case of 
westerly winds, St. Louis River water Fig. 13), was transported along the 


‘needed ‘to drive contaminants away from the ‘south shore and trap them in 
the southwest corner of the lake, 
— Direct verification of the long-range dispersion of runoff cannot be made, 
_ since it would be cost prohibitive to run the model for thirty days of actual — 
OM winds until the next useful Landsat image is available. However, an indication — 
ris Of 5 whether the modeled results in Figs. 12 and 17 are realistic can be made a 
consideration of Landsat data. Using criteria developed from analysis 
of optical properties of the contaminants and correlation of Landsat multispectral — 
a data with ground truth observations (17,18,20), the distribution of runoff three — 


weeks after the 1975 peak runoff of the Nemadji was derived (Fig. 14). The 


FIG. | Louis River Effluent and Red Clay Concentration after 
1975 Spring Runoff Obtained Using Landsat Data for May 18,1975 re 
1975 runoff occurred later in the se season, and was larger than the 1976 runoff. y 


However, the winds following the 1975 runoff were variable, and produced 
a distribution of red clay and St. Louis River runoff quite similar to the simulated 
distribution of rayon contaminants shown in Fig. 12 for comparable time and > 

A rough anes on — validity of the ‘omman runoff dispersion after one _ 
- month can be made from examination of Landsat data for early May 1976. 
Fig. 13 shows a comparison of the simulated dispersion with the actual plume — 
distribution on May 3 1976. The remnant red clay plume due to runoff, indicated _ 
in Fig. 13 by the turbidity hook (solid line) north of Brule River, is quite comparable _ 
to the simulated runoff plume shown by the dotted broken line. Both hooks ; 
have roughly 3 mg/L concentration . of red clay. ‘Comparison of 
in Fig. 13 near Duluth is not meaningful because strong northeast winds in 
late April ihe generated significant shore erosion and een The simulated © 
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plume density and the actual plume de density { for May 3, 1976 differ north of 
the Amnicon River because the simulated plume accounts only for turbidity © 


west of Brule River is still quite > evident in the May ry 

_ data and simulated results show that the lower contaminant concentration contours 
are confined to the area of the lake west of the line between Knife River, 

_ Minnesota, and Brule River, Wisconsin (Fig. 1). Comparison of lower concentra- 
tion is, however, difficult because of the ‘uncertainty in effluent identification. — 


The Lagrangian marker solution technique worked well in simulating both 
the development of the Nemadji red clay turbidity plume and the long- range 
4 ee of runoff water from the Nemadji and St. Louis Rivers. This technique Z 
_ yielded bester results than the straightforward application of mixing and diffusion 


in finite difference methods. Currents from a depth integrated transport model 
_ with a relatively coarse 3.7 mile x 3.7 mile (6 km x 6 km) grid yielded accurate 


with a finer grid in western Lake Superior would further improve the accuracy 
_ Numerical simulation of the transport of pollutants also answered many 
Gas which arose in monitoring the 1976 spring runoff. The origin of many ; 
3 _ distinguishable features of the Nemadji red clay turbidity plume became clear. | 
The relative contributions of Nemadji runoff, resuspension, shore erosion, and 
| runoff from south shore streams could be estimated. We estimate about 65% © 
of the red clay turbidity observed in the Landsat images of April 5 and April 
. 6 1976 came from the Nemadji River. The distribution of contaminated water _ 7 
_ in Lake Superior for various wind conditions was also explained. For westerly © 
winds, effluents from the Nemadji and St. Louis Rivers are carried out of — 
western Lake Superior along the south ‘shore. ‘The “onset of easterly winds 
_ generated eddy circulations which tend to trap these river r effluents i in the western _ 
arm of the lake. The model indicates a conservative contaminant from the St. 
_ Louis and Nemadji Rivers might not be flushed from the Duluth | —_ until — 
_ the se severe storms of November amber change the the ‘prevailing « current patterns. a ato 
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plume simulation. This accuracy in turn confirmed the reliability of the current q 
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The following symbols are used in this paper; der 
= size of marker particle before diffusive change; 


depth averaged concentration; 
= background concentration for individual marker particle; 


= 


H = 
= of surface of water; 


Zz - 
_R* = first approximation of new marker particle net transport position; q Yo 
marker particle position after translation by current; 


averaged x component of f current; 

depth averaged y component of current; 
‘net transport current at Rows 


net transport current at R*; 


v= = three-dimensional gradient operation; —__ 
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DISTRIBUTION OF CREST-TO- ‘TROUGH WAVE HEIGHTS 
By M. Aaiz Tayfun," J ASCE 


—_ It is known that the dintribution of wave heights in a sea state is — 
; _ well by the Rayleigh law under certain conditions. Specifically, it is assumed 
= the sea surface is linear Gaussian ' witha narrow-band spectrum. However, 
it is also known that the characteristics of a realistic wave field are not always 
consistent with these assumptions, and various discrepancies are often noted 
: between field observations and the Rayleigh theory. Likely causes for these 
_ have been considered in a series of studies by Forristall (2), Nolte and Hsu 
(8), Longuet-Higgins (6), and the writer (9, 10). They include: (1) The nonlinear, 
~ non- Gaussian characteristics of the sea surface; (2) effects associated with 
_ wide-band spectra; and (3) wave- breaking in deep and shallow water depths. _ 
_ A further cause for the discrepancy between the Rayleigh law and empirical 
_ data is related to a difference in definitions. The Rayleigh distribution is especially — 
relevant to wave envelopes, which can be viewed in time or space as a pair 
of symmetrical curves passing gi wave crests ome troughs, 


so 0 that maxima and minima of the upper and lower envelopes exactly éoincide 
with wave crests and troughs. At the other extreme, if the sea state is characterized 
with a wide-band spectrum, the wave profile and its envelopes vary rather 
\@ - _ rapidly, and local maxima and minima of the upper and lower envelopes do 
necessarily correspond | to wave crests and troughs [see, e. Ref. 
either case, the difference between the upper and lower envelopes, or, 
equivalently, twice the local wave amplitude, represents a “wave height which — 


g Almost all empirical ‘distributions derived from field data are based on a 


crest-to-trough representation which differs fundamentally from the envelope- 
_ based definition. In fact, a crest-to-trough wave height corresponds to the sum 
_ of two values of the upper envelope s separated | by approximately half the local ‘ 
i period or wavelength, and represents an averaging operation on the envelope. 
Therefore, the objective here is to examine explicitly the distribution of pain 
heighes based on this latter definition and to compare it with the Rayleigh law. 


_ ‘Assoc. Prof., Civil Engrg. Dept., Kuwait Univ., P.O. Box 5969, Kuwait, 
4 Note.—Discussion open until January 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on May 6, 1980. 
This paper is part of the Journal of the Waterway, Port, Coastal and Ocean Division, 
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_ we aim to show that the density of crest-to-trough wave heights is fundamentally _ 
different from the Rayleigh form, showing 2 an \ excess of waves with heights — 


= larger values of the spectrum bandwidth. 


or Wave Heicuts Basep on Wave Envetore 


z the elevation of the sea surface at a given point be represented by — 


€, £08 (0, t+e 


= uniform 


random over an interval of 27; and, c, such that over 


S(w) d 


with S(w) representing - one- -sided one- -dimensional spectrum. The jth spectral 


In particular, Lo = and w, = = mean The 

spectrum is considered to be narrow-band if (4,5) 


In general, n(t) can be re rewritten as 

= cos (wot + 8) (5) 


in which . A (t) and = = the functions, ‘Tespectively, 


exp (10) = exp t+e, +6 


in which 4 =) sin (w,t+e€,) 
The free surfaces, and 4, become asymptotically Gaussian and independent 
s the number of terms, N, of cosine functions approaches infinity. ‘Under 
this limiting condition, the probability density of A is the Rayleigh law (10). 7 
4 In particular, the density corresponding to = rece amplitude, € = A/A,.. 
The joint density of = E(t + is Isee, ¢ Ref. 
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in which ie o(* ) = zero- ee a of the first kind; and 
with = 


Bove 
wd lower wave envelopes (9,11). The basic character of these is intimately © 
related to the spectral bandwidth measure, v. In sc al = dA/dt is z 
zero-mean Gaussian with variance, v i.e. 
| Thus, the narrow-band condition v? << 1 suggests that the underlying spectrum, J 
‘S(w), and f ,(x) are centered sharply around w = w, and x = 0, respectively, 
and behave as pseudo-delta- functions. Correspondingly, each realization of y_ 
ma depicted by Eq. 5 has the distinctive form of slowly modulated | oscillations | 


crest, 7 =0 and AU = . Similarly, when 7} is maximum, corresponding to _ 
a trough, y = 0 and A = %. Hence, we have the upper and lower envelopes, 
ei —A, from a pair of symmetrical curves which pass through wave crests” 
and tro ns espectively, and va lowl between as implied by E 
n sughs respe tively, an nd vary very slowly in between y 


of phase by sc so that when is ‘maximum, ‘corresponding to wave = 


J _ When the spectrum i is with rae so that the condition v” << 1 is not satisfied | 


_ Strictly, Eq. 14 suggests that the rate of change of A could be rather | large. 
In this case, realizations of y can no longer be viewed as amplitude- modulated ; 
oscillations with an apparent mean frequency, w,. Consequently, the physical _ 
interpretation of A as an is somewhat lost. Nonetheless, the 

stil be defined as = “2A. Noting that H,, it then follows that 
the scaled height — = -H / Hm. has the invariant probability density depicted 
cx _ by the Rayleigh form of Eq. 9, Tegardless of whether the condition v? << 


1 is satisfied or not. 


ok 


be The cre crest 10. rough wave height i is defined 


in which A A(t) and A, 2 = A + 7/2) respectively, two values" 
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of the w wave by half the lo local period to 


Fal or 
e conditional distribution T= = T/2, is defined by 


te 


Finally, the period tot take all possible values, we 
probability density of 


which f(t) = the density of 7, with ¢ = 0. 

Since the joint density is closely related to the 
et 10-13, the probability “density of crest-to-trough heights is a a functional — 
of both the spectral shape, S(w), and the density of wave periods, ta This 
i essential characteristic is a distinct contrast to the Rayleigh law. It also renders 

_ the explicit evaluation of f; numerically a rather tedious attempt in ae 
_ Therefore, i in n order | to proceed with this Teasonable ease, Eq. 


<< Under this condition S(w) as a pseudo- function 
centered around w = wy, and the Tayler series 


be used in Eqs. 12 and 13 to thet = 


| 


| 
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4 _ Making use of Eq. 10 and paying appropriate attention to the region of integration = &g 

implied in this definition, we can writeitas 

| 

} 


Proceeding further, the corresponding density, f (t), is ety centered around — 
the mean apparent period, 7, = 27/w,. Therefore, it behaves as a ceiiiediine. 


_ function with respect to the integrand in Eq. 20, reducing | Rfesisaboinl aff ; 


= 2 —u,u; - du . 


od: ne woile ba snes 


om vad 
“The quality of this approximation should be rity satisfactory provided 
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FIG. ahaa between F Rayleigh Law and Density oll 
Heights under Narrow-Band Conditions 


‘picture of the basic behavior of f;, which still has to be determined numerically 


2 but with far more ease and less information than that required by Eq. 20 m. 


_ Using Simpson’ s rule, ‘Eqs. 10, 23, and was evaluated 


1, corresponding to v = 0.1, r = 
0. 95, and v= 0. 2, = 0. 80, respectively. The Rayleigh density is shown 
inthe same f figure he comparison. Iti is observed that the density of crest- eee all 
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5 heights s displays an excess of waves with heights near near the m midrange an anda Sittelaties: —_ 

j 1 at the two extremes in a manner increasing with v. As previously suggested - 

in the Introduction, this is because the crest-to-trough height as defined di 
Eq. 16 represents an averaging or smoothing opesetion on the wave envelope, — ae 
damping its oscillations about an average value. Thus, the resulting density 4 4 
must have the s same mean, and show an excess about the mean and deficiencies — a 
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FIG. of Scaled Crest-1 to- Wave Heights from Bretschneider ( 

_ in Comparison with Rayleigh Law and Theoretical Density of Crest-to-Trough Wave 

_ The theoretical mean and variance of € can be obtained exactly, by utilizing — 
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Since r? always | (except when is monochromatic wave with a constant. 
= it follows from a comparison of Eqs. 26 and 27 ‘that o 
which is consistent with the preceding analysis and Fig. 1. 
- _ The ultimate effect of a deficiency towards the high wave tailis an ‘underestimate 
] in the exceedance probability of larger waves. The partial reason for this kind 7 
of discrepancy systematically observed between the Rayleigh form and field — 
_ data can be explained on this basis. A case of practical interest has been considered 
7 an by Longuet-Higgins (5), based on fieid data collected in the Gulf | 
of Mexico. The data was compiled earlier by Bretschneider (1) and contains 
_ approximately 400 crest-to-trough wave heights in records e-11-3-15 of the same _ 
__ ‘Teference. The histogram of crest-to-trough heights in a scaled form was | 
reconstructed here in Fig. 2, using the scatter diagram given in Fig. Sivce 
Ref. (1). The theoretical crest-to- “trough density of Fig. i with v = 0.2 and > 
the Rayleigh form are shown also in Fig. 2 for comparison. It is evident that ' 
_ the histogram does display the discrepancies under discussion relative to the is 
Rayleigh form. This was noted also by Longuet-Higgins (5). On the other hand, 
most of the histogram ordinates are predicted better by the theoretical zero-cross- 
ing wave-height density corresponding to v = 0. 2 which closely approximates — 


the 0.234 estimated by Longuet-Higgins (5) for the same data. Discrepancies — 


of this sort can be seen in other empirical results such as those given in Figs. 


/4,21-23, 4.28, and 4.30 of Ref. + 
Conc.upinc 
_ The definition of crest-to-trough wave heights differs from the envelene “based 
- representation relevant to the Rayleigh law. It was demonstrated that, 
J comparison with the Rayleigh form, the density of crest-to-trough heights pra Aer 
an excess of waves, with heights about the mean wave height and deficiencies _ 
at both ends away from the mean. A principal effect of a deficiency towards __ 
the high wave tail is an underestimate in the exceedance probability of larger 
waves associated with the greatest engineering concern. Considering that the 
representation of wave heights based on the wave envelope is consistent with 
the theoretical and physical specifications of narrow-band sea waves, there is 
“no obvious justification for using a crest-to-trough definition. This plea is 
_ Supported ultimately by the need to explain the differences observed between __ 
field data and the Rayleigh theory. The rationale is that, if future field data 
is also analyzed on the basis of an envelope definition, possible discrepancies 
arising from the crest-to-trough definition would be eliminated, and, thereby, — 


effects Physical mechanisms such as wave breaking could perhaps 
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100-150 sample points for the construction a 
_ In principle, a sample of this size is insufficient to estimate a stable, well-resolved 7 
distribution. ‘Very often gaps exist in the ‘data so that the ; analyst is forced 
- to enlarge intervals. The histogram displayed in Fig. 2 would actually have 
no ordinates in several intervals if the mesh presently used for the abscissa 
_ were halved (see, Ref. 1, Fig. 5.1). For the same record, the cavelegs presents 
a continuous sample space with a much larger size. To be specific, a 20-min 
record sampled at a periodic interval of A = 0.5 sec would yield 2,400 | sample 
- points, which would allow more flexibility in performing data analysis and 
obtaining histograms ordinates with varying degrees of statistical quality. _ A - 
_ The construction of sample points from an envelope is simple, and can be _ 
“carried out efficiently if coupled with FFT operations, routinely used in the 
ent analysis | (of sea surface elevations. Consider, for instance, the finite — 
record Ths in which n, mtd). The FFT coefficients are given 


exp 
bp 


= complex conjugate gate of Now define OF.» 


H H, = 2(n? + 2). = 0, 1, dive 


= crest- to-trough wave heights from different records are grouped —— 


for Rayleigh- -distributed variables virtue of the invariant form of this distribu- 
tion if the variates are scaled with respect to their means or root-mean- square. 
- values. However, if samples are not from a common probability law, it is obvious 
that the histogram constructed on this basis cannot be entirely meaningful. 
Therefore, | any other physical reason associated with a sea state put aside, 
the reliance on the sole usage of crest-to-trough wave heights inherently requires 
that samples from different records with varying spectral bandwidth measures 
not be grouped together since they do not ay the same probability law. peat x 
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The proordure (EX Wi") prevented here allows for wave deformation 

The following symbols are used in this paper: how ¢ the design 
flesitaiity m ‘the chuive of tbe Gove, wheres: the abiliqnic Weat 


A = wave envelope or amplitude; of wave 
A,,A, = values of A separated by half wave period; 


= Fourier coefficient of nth component wave; _ iy tai 


= expected value operator; in the frequancy aad 
= conditional | density of é, given T; vehagity poten thet 
F, = hypergeometric function; of the Roundary 
. probability density of wave period; 
conditional density of given 7; 
density of and &, respectively; 

wave height based on wave Gn 

crest-to-trough wave height; Const and 


\d discrete values of H and H, respectively; "107, 


spectral density; 
wave period; ant 


| 
| 
| 

a 
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= integration variable, in Eqs. 19, 20, and 24; hate 
random phase associated with the nth component wavelet; = 
sea surface elevations as a function of time; 
= discrete values of and 7 respectively; 


= = nondimensional crest-to- -trough w wave height; 


* = values of separated by half wave ‘Period; 
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_ EXTENDED VELOCITY POTENTIAL WAVE KINEMATICS 


By Kostas F. Lambrakos' T 6,1} ime (270) 
S521, odie bemoqe: sqade av avew Vo 
_ Wave theories are an integral part of many design nimeeieie for offshore 
is and coastal structures. They are used with the design wave profile to predict 
flow quantities such as wave velocity, acceleration, and pressure from which 2 
design force calculations are made. _ Wave theories have also played a role in df q 
the determination of the e empirical coefficients in relationships that translate 
wave velocity and acceleration to wave force, e.g., Morison’s equation. In such 
_ applications, measurements of force are correlated with wave velocities and = 
- accelerations predicted by wave theories from measurements of wave profiles. 
Several nonlinear wave theories (2,3, 10) have worked the problem of symmetric 
harmonic waves (Stokes they have been used extensively by engineers 


‘maaied for measured waves on the basis of the zero-crossing wave height x 
and period. The stream function procedure (3) can treat observed nonsymmetric ; 
= waves by achieving a best fit to the measured wave profile with — 
a Stokes-type wave profile shifted i in time (or space). The key assumption made # 
in these procedures is that ocean or design waves propagate unaltered in form 
(shape). As a result, their domain of applicability includes only single waves _ 
of fixed wave-crest-to-wave-height ratios for a specified wave height, period, a 
_ The procedure (EXVP) presented here allows for wave deformation (a charac- 
teristic feature of ocean waves) and its domain of applicability includes, in ; 
7 ‘ addition to Stokes waves, , waves of various crest-to-height ratios and groups _ 
: of waves. The ability to vary the crest-to-height ratio allows the design engineer _ 
_ ‘greater flexibility in the choice of the design wave, whereas the ability to treat 
irregular waves may be used to improve the prediction of design wave forces © 
through better correlations of measured waves and wave forces. 
‘The EXVP procedure solves the wave velocity problem with a velocity potential 
function defined he a double Fourier series expansion in the frequency and 
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(OTS) experiment (1,4,6). The procedure has also been used to | investigate the 
dependence of wave force on wave shape as reported at the Offshore 
Conference (9). Additional theoretical analysis on EXVP ‘and 


—The water motion waves at sea is 
defined on the basis of Fig. 1. The irregular wave profile extends indefinitely 
in the x direction, in water depth, d, over an impenetrable seafloor. With the | 

FIG. 1.—Definitions for Two-Dimensional Wave Flow Problem of 
assumptions of two dimensions, irrotational flow, and incompressible water, 
equations of motion for the wave can be in terms of the 


continuity: — —s 


z+— + (uw? + ——— = (2) 


in which x, z = the horizontal and vertical coordinates, respectively; p = the 
by seca pressure; p= = the water density; g = the gravitational constant; a 
and c = an undetermined constant. The horizontal and vertical components, i 
u and v, of the water particle motion are related to the velocity potential function, | 
by its conve derivatives as a consequence of the irrotational flow assumption: 


/ 
fitted numerically. The intrinsic assumption § aa 
ure is that waves are periodic in time and 
space—an assumption which is also made by all Stokes-type procedures. 
ge _ This paper describes the theoretical and computational features of the proce- _ — 
| 
- q 
: 


sea floor over ali all space and time. There are _ boundary conditions 
wave surface and one at the sea bottom: 


m th Ww d d b 


The dynamic boundary condition is Bernoulli’ s s equation evalu ted 


which p, = the atmospheric pressure. dod vd vit 
The boundary condition at floor that the vertical | component 
of eee be zero, 


In view of Eq. 3, the solution to the wave velocity problem previously ‘eed 


consists of finding a velocity potential ' which satisfies Eqs. | and 2 and boundary ~ 

7 conditions, Egs. 4, 5, and 6. There is difficulty in - obtaining such a 1 velocity — 

in potential since the boundary conditions at the wave surface are nonlinear. 
Consequently, it becomes virtually impossible to satisfy them exactly, pase a 
when the wave profile is not known a a priori over space and time, as is — 


which th the ‘commonly known wave ‘theories make, reduce the diffi iculty of 
problem as previously stated since the boundary conditions need only be satisfied 

for one space point and over the period of the wom, Since profile deformation — ; 7 
is a characteristic feature of almost all waves, the wave velocity procedure us 
presented here attempts to satisfy the hydrodynamic equations and boundary — 
conditions without the restriction of constancy in the wave profile. uo bsivtu od 
_ Extended Velocity Potential (EXVP) Solution. —The ‘velocity potential of the — 7 
present approach is assumed to be given by a Fourier Series expansion in the — 


z,t)= cosh (k, cos (Kk, x) x)- + F, sin (k,, 
in which | an F, = functions of time; k, = the wave number of the nth Fourier | = 
component; and E,, F,, and k, are evaluated from the boundary condone 


satisfies exactly the pete esta equation and the boundary condition at the bottom. a 
i, It is believed that this general expression could be used to obtain a numerical 


| 4 
= soluuon [for the case where the wave promue iS specirie a ore han One = 


s i and F, are defined appropriately, the velocity potential — 


terms of progressive harmonics: ane 


in = constants. This: form of the | potential function 

“differs from previous definitions in a crucial way. For each frequency, w,, 

. ‘ea is a corresponding set of wavelengths, {k,}, and vice versa. Consequently, ; 
~ the harmonic wave components for a specific frequency value are allowed to 
propagate with different phase velocities defined by the ratio, w w, Jk,. This feature © 
a EXVP enables treatment of nondeforming | (Stokes) as well as ; deforming 
wave shapes. The proposed method differs from other procedures that treat 
the nondeforming wave shapes based on one-to-one correspondence between 
frequency and wavelength, and on the same peypagy to-wavelength ratio 

constants, k,,, By», in the EXVP expansion, Eq. 8, are determined 

_ numerically by satisfying both boundary conditions, Eqs. 4 and 5, at the wave 

_ surface. Within the assumption of wave periodicity in space and time and with 
“the specified pressure at the wave surface, the harmonic coefficients, A,,,,, B rm: 

that are dominant in the expansion depend on the wave profile. If the wave 

profile is the Stokes shape, the dominant harmonic coefficients are only 
those with equal w wave frequency- -to-wave number ratios, and their values are 
the same with those for the Stokes-type solutions (10). If the wave profile 
differs from the Stokes shape, additional terms become important. In the case — 
of a relatively small wave, the coefficients corresponding to Airy wave fre- 
quency-to-wave number ratios dominate. 
The EXVP ‘Computational Pri Procedure.—The constants in EXVP are e calculated 
by an) iterative numerical procedure that minimizes the error in the kinematic 
and dynamic boundary conditions. The specified wave profile, which can be 
of more than one wave, is assumed periodic in time and space. Thus, the 
numerical fit to the boundary conditions is for an appropriate (and finite) number — 
of points in space and time. Since the wave is deforming, matching to the 
boundary conditions is carried out over both distance and time; if the wave ; 

_is nondeforming (same celerity for all Fourier components), the matching may — 
carried out over only distance or time. | 

_ Since the wave profile is normally specified at one space location, the error 
in the dynamic boundary condition (Bernoulli’s equation) is computed for only © 

that one location. For every | other location the error is zero since Bernoulli’s 

_ equation, Eq. 2, is used to calculate the wave profile for any location other 

than that for which the profile is known. At the end of the computation, a 
_ new wave profile may be calculated from Bernoulli’s equation even for the — 

- location at which the profile was given; this wave profile should approximate 
well the initial profile and, along with the velocity potential coefficients, would 
AS was stated earlier, the EXVP calculation a assumes es that the specified wave 
profile is periodic in space and time. The fundamental period i in time is normally 


equal to the time span, 7,, of the specified wave profile. ~The frequencies, 


| 

can be express — 

— 
— 

(x, z,t) = 
> 

q 

{ 

| 
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‘tk for the Fourier series expansion are defined by b: 
in which w, =2n/ T, . The unknown fundamental period, L ,,in space is calculated _ 
by” the error minimization yn procedure; the corresponding wave number is ky “4 
= 27/L,, and the remaining wave numbers are defined by 


The unknown coafficients, (B, in the EXVP expansion are matrices 


each element of which to specific values of and wave 


er i.e. — 
{A am} wine = An gigs de? 


in n which the rows correspond to the variation in the frequency, ‘subscript m, 


and the columns to the variation in the wave number, subscript m, and each - 
frequency and wave > number is given by Eqs. 9 and It 10. The diagonal elements 


in cases of nondeforming waves. Even for deforming profiles, not all eniaiae 
_ in the coefficient matrices need be evaluated for a satisfactory fit to the boundary 
conditions, , and the judicious choice of elements can save in computer time. 
The constant, c, in Bernoulli’s equation, Eq. 2, is also considered an unknown. Eo 
- The: numerical procedure calculates the values of the unknowns in the = : 


@ expansion ‘that n minimize the error in the boundary | conditions in a least-squares 
4 sense. With given values for the unknowns, errors in the kinematic boundary — 
_ condition are calculated for every point in space, x,, and time, ¢,; for the dynamic 
_ boundary condition only errors over time are calculated for the point in space 2 
for which the wave has been specified. The points may be given” 4 


Ax and At = half the shortest | wave and ¥ wave > period expected in the 


The errors in the A conditions are calculated as s follows. the 
in which the right- mane ‘sides of Eqs. 12 and 13 are eva evaluated at the wave 
7 surface for time, ¢,, and location, x,. The least-squares sum of these errors 


q 
: 
— 


This is the quantity which is minimized to obtain the ¢ optimum values ‘for the 
ae the subscript, p, designates the spatial location for which the wave 
profile is given. The minimization is achieved with a variant of the Gauss-Newton 
_ technique (5), which iterates on the first estimates of {A,,,}, {B,,,}, c, and ra 
S ,}. A measure of quality of fit to the boundary conditions is obtained by _ 
(weighting) the errors ,(E, and by the root-mean- square 
values of the wave profile and its time ‘derivative, respectively 


The first: estimates of the coefficients are arbitrary but judiciously ‘picked, 
* will be demonstrated in the examples given as follows. Some of them may © 
“f be arrived at on the basis of linear (Airy) wave theory. Depending on how > 


close the first estimates are to a good solution, the computer time | ee = a 
if instead of the wave ‘Profile, the wave height and period of the wave are 
“specified (for a given water depth), a variation of EXVP (EXVPD) | can be 
used to calculate the wave shape. This case is the most common in the design 
_ of offshore structures, and has been treated with the Stokes procedures. For _ 
the EXVPD solution it is required that only the diagonal terms in the coefficient _ 
= matrices be retained and that a constraint ¢ condition in the wave height be included | 
as part of the least- -squares ‘minimization procedure. The initial estimate of the 
ue shape can be simple harmonic or a wave shape derived from linear or y 
Stokes theory. In the end, EXVPD provides a wave that propagates without — 
deformation and has the specified wave height. 
‘The time requirements for executing EXVP in the IBM 168 computer wre 
vary from several seconds, for | the nondeforming waves, to several minutes 
for a wave segment of two ‘irregular waves. In the final analysis, the time 
required in each case depends on the number of coefficients used, on how 
_ compatible these coefficients are with the given wave profile, and on the number 
of points in space and time required for fitting the boundary conditions. With 
present computer capabilities and EXWP’s current numerical procedure, the 
aa 
_ For the purpose of illustrating some of the features of the EXVP solution 
and how it compares with field measurements of waves and velocities as well 
as with the Stokes type procedures, the EXVP procedure has been a : 
to two simple wave profiles and several irregular ocean wave profiles from __ 
the Ocean Test Structure measurements. The two simple wave profiles were — 
ic of identical wave height, 40 ft (12 m), and wave period (12 sec), and the waier — ; 
depth was 100 ft (30° m). One was a cosine function of time and the other 
a Stokes type wave profile. The irregular wave profiles were composed 
= of two adjacent waves, one with a dominant crest; their crests and troughs — 
are given in Table 3. The water depth was 66 ft (20 m). | ei Ose “= 
Cosine Wave Profile.—The initial and final estimates of the EXVP coefficients _ 
for the are given in in Table 1. The improvement in the initial estimates 


i 
| 
Ne 


bis suggests that for a good solution the first guess 
_ does not need to oe : very close to the final solution. It is also evident from _ 
§ _ the final estimates that there is significant interaction between the various wave _ 
& Components | since the off-diagonal terms are not “negligible. The interaction 
manifests itself in the deformation of the wave profile and also in the sane 
- velocities (8). During propagation the cosine wave deforms to wave forms “a 
_ significantly higher crest-to-height ratios than the ratio of the Stokes wave of 
_ the same height and period. The near-surface water particle velocities reflect 
_ differences in the e wave crest—the velocities for the steeper waves are” 
higher than the velocities for the Stokes wave. 
_— Stokes Wave Profile. —The EXVP solution converges to ‘the Stokes velocity 
a potential solution for the same as Stokes waves. The 


-300 | 10:84 


Stokes wave me calculated with EXVPD for a wave re period, and 
water depth of 40 ft t (12 m), ‘12 sec, and 100 ft (30 m), respectively, was used 
ith EXVP, and the results are shown in Table 2. The final EXVP coefficient 
matrices in Table 2 are dominated by the diagonal terms in the {B} matrix, 


ocean test structure experiment is very similar to the treatment of single waves. . 

The OTS measured simultaneously storm wave profiles at and water particle 
- velocities at four elevations in the water column and, consequently, provides — 
a unique opportunity for testing the velocity predictions from EXVP. soap. ol 
is the OTS data include a great number of waves, only 10 selected irregular waves — 
Boxee were used for this study. These waves, however, are sufficient to demonstrate 
_ some of the EXVP af the EXVP capebilities onl sea waves and to compare these capabilities 
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on 
are very Close to the velocity potential COeIlicients OF the STOKES SOLUTION. 
aw, Pra ac The P treatment of irreoular waves am the 


_ The wave profiles i in this study were composed o of two adjictat waves, whose 
troughs, crests, and time span are given in Table 3. No particular \ criteria were 
used to choose these waves other than that the waves be large and cover a % 
range of wave heights and periods. However, some bias toward two-dimensional 
2.—EXVP Matrix Coefficients for Stokes Wave 


INITIAL ESTIMATES FINAL ESTIMATE 


pe Note: Stokes velocity potential coefficients: 877.50, —122.30,8.84. = = © 

TABLE 3.—Irregular Ocean Test Structure Waves = = 


WAVE CREST 2 WAVE TROUGH ‘wave TROUGH 2) WAVE CREST 1 | WAVE PERIOD ’ 
(IN FEET) 
= 


= 


NOTE: 1 FOOT = 0.305 METERS 


6 
-10.00 | 00 | 00 | 00 | -3.72 | 00 | 00 
bey 5.00 4 | 00 | 00 | 261 149 | 00 00 
--3.00 | -3.00 | -3.00 | 00 | -040 | 3.08 |-016 | 00 = | 
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-3.00 | -3.00 | -3.00° 70.01 vy 4.30 | 8.75 
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0.0 3.00, | -3.00 3.83 | -0.85 | 
a 0.0 0.0 | 2.00 | 200 0.0 0.0 | 1.33 
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the swave crest of interest. The study focused on the velocities due to the = 
tanger of the two zero-crossing waves in the wave profile and compared the 
extreme (crest and trough) values of the measured velocity magnitudes for the 3 
S crest and preceding and following troughs with the corresponding peak values 
predicted from the EXVP and Stokes (EXVPD) procedures. For one of the : 
_ waves, no. 10, the study aiso compared the entire time record of the wave _ z 
profile and the velocity magnitude for every elevation. Although the velocity 
4 direction changed within each wave prpent, no attempt was made to take direction - 


_wave SEGMENT___, 
"FITTED WITH EXVP 


= 


\ 30 


Le 
.~@y 


CALCULATED (EXVP) WAVE SURFACE 


oy 
—EXVP Wave Profile Fit to oTs- Measured | (Wave No. 10) 


a account since the extreme velocities, which are of special interest, were 


nearly collinear. The velocities, however, were adjusted for the component 

of steady current (hourly velocity averages) which was present in the data. 
Also, the extreme velocities chosen for comparison were those that occurred _ 

- within +0.5 sec from the instant of the maximum (crest) 0 or the minimum (trough) - 

; a While in the EXVP analysis the fit to the wave profile was carried out in 

= terms of the actual wave trace, the EXVPD analysis was performed on the 


basis of the zero up-crossing and zero down- ~crossing wave heights and and 
% 


4 | 
— 
The 
5, and 
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which are the c common approaches of approximating wavesat sea awith Stokes-type 7 


“used fe for the simple wave calculations above; it was assumed that the four 

frequencies would be sufficient to describe the wave velocities. 
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FIG. 4.—EXVP Wave Velocity Magnitude as to OTS-Measured Velocity 
for Various Elevations (Wave No. 10) sone 
< The comparison on of the EXVP solution with the entire measured wave and _ 
velocity traces has been carried out with wave profile no. 10, which has 
- substantially different troughs, 6.66 ft (2.03 m) and 4.08 ft (1.24 m). The 30-sec a 
wave trace from which wave profile 10 was chosen is shown on Fig. 2. ‘Fig. 

; 3 shows the two-wave segment of this trace that was used in EXVP and the — 
& that represents ths EXVP fit to this segment. The EXVP profile closely 
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due to the lack of higher | higher besnidtend in the velocity qutattieh: These differences, 
however, do not strongly affect the velocities since the measured and calculated _ 
_ velocity profiles for all lower three elevations are in very close agreement, 
as shown in Fig. 4. For the lowest elevation, a small phase shift between 
measurements and calculations for the crest velocities is present, but its origin 
has not been determined. For the elevation above the mean water level the 
_ measured velocities are not very reliable, due to poor current meter Tesponse — 
7 as the sharp spikes and * “‘ringing’’ in the data indicate. However, for this elevation 
_ the velocity predicted is in generally good agreement with the measured velocity Pe } 


_ For the same elevation, the calculated velocity for the first crest in the wave 

: profile is peer since the calculated crest elevation was less than the highest 7 
onl The Stok es wave approximations to ‘the wave of interest are shown in Fig. 

4 5, and the measured and Stokes velocities for the various elevations are compared © 
in Fig. 6. It is evident from Fig. 5 that the zero down- crossing =e matches ill 
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3 FIG. 5.—OTS Wave Profile as Compared to to Zero Up- Crossing ai and ea “ iii 
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owing ‘to the ur uncertain quality ‘of the “measurements, no definitive comparisons 
can be drawn. The calculated velocities at the 51-ft (15.6-m) elevation from 


trough, by as much as 25% and 10%, For the lower 
a the calculated velocities a better with the data, the up-crossing velocities 
ies, and closer to the data. 


+ 
. 
7 
a profile, whereas the zero 
-up-crossing wave differs significantly from the corresponding measured wave §& 
- profil in both the crest and the trough. However, the calculated velocities 
_ for the two Stokes waves show the opposite trends, i.e., the up-crossing velocities 
- | are closer to the data than the down-crossing velocities. The calculated velocities _ 
hath Stakec waves are overall larcer than the data for the crest and following 
_ The crest and trough values for the zero down-crossing wave are larger 
: - those for the zero up-crossing wave, and the same is true for the velocities. _ 


i noted previously, ‘the up-crossing velocities ore in better agreement w 

3 the data for these particular wave profiles, but, for other irregular wave profiles 
for which the down-crossing wave is smaller than the up-crossing wave, it is a 
conceivable that the down-crossing velocities will be closer to the data. The © 

_ differences between the velocity measurements and the Stokes velocities are ; 
greatest for the following trough for which the eerrpansrs between celoaiated 


measured elevations | are 


ELEVATION = n hes aa 

MEASURED a 

=== ZERO DOWN-CROSSING 


D) 
~ 


on 


VELOCITY (FEET PER SECON 


4 ELEVATION = 31 FEET 


ELEVATION = 16 FEET 


ay 10 12-14 16 18 20 
(SEC) 
‘wie’ 1 FOOT PER SECOND = 0.305 METERS PER SECOND 


FIG. 6. —Stokes Wave Velocity Magnitude as Compared to Velocity 
for Various Elevations (Wave No. 10) 
e Stokes wave would match both the crest and trough elevations of the measured | 
wave profile, the Stokes velocities would, in general, be larger than the ieee | 
_ velocities. These observations suggest that the accurate prediction of velocities 
q for irregular waves using Stokes representations is not a straightforward exercise; - Ks r 


even if the ‘‘right’’ height is chosen for the crest velocity, the trough velocities 


the of the EXVP and Stokes procedures, it becomes evident 


that the EXVP solution matches the data closer than the Stokes solution. This — 

is true for the extreme crest an and trough values as well as the rest of the wave = 

and velocity p profiles. The accuracy of the velocity ‘prediction before and after 
the wave crest can have a twofold impact on the force calculation since the 
maximum force, e.g., for a vertical pile, normally occurs away form the wave ~ 
crest, and, thus, both the flow velocity and resulting acceleration contribute 
i. the force. It should also be noted that the Stokes procedure fails to ea 
higher velocities under the preceding wave trough than under the wave crest, _ 
as is observed in the data and predicted by EXVP. 
_ The measured extreme velocities for the 10 waves studied have been differen- 


tiated on the basis of the crest, preceding trough, and following trough values 
_ and in terms of the four elevations in the water column. These data have been 
: compared with the corresponding data from the EXVP and Stokes representations _ 
(Figs. 7-9). be. the 10 0 waves only nine waves were used in these figures 
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7. —Comparison of OTS- Measured Extreme Horizontal 


since the extreme velocity values of wave no. 2 were ae the time band» 

+0.5 sec about the corresponding extreme values of the wave profile. The — 

measured values for the highest elevation that were used in the comparisons: / 

e only those that were subjectively judged | of good quality, = Saree 

_ Overall, the -EXVP velocities, Fig. 7, are in better agreement with the 

measurements than the Stokes velocities, Figs. 8, 9; this is especially so under 

the wave troughs. The EXVP correlations show no significant differences in , 
¥ _ terms of elevation; the crest velocities show minimal bias (on the high side) 
a and about 10% scatter; the velocities for the preceding trough also show 1 minimal 


> if any, but more scatter, about 20%; the velocities for the following trough 
show a small bias, on the low side, and greater scatter than the velocities 
for the preceding trough. Part of the scatter present is probably due to wave 
spreading. However, the differences between the crest and trough velocity 
correlations are consistent t with» the differences between the measured and 


7 


‘for the troughs are ums than those for the crests. It is s believed that if the 7 


| 8 


me 
matching of the trough elevations were to be improved and a greater ‘segment 
of the measured profile were taken into account, the scatter in the trough velocity ef a 
correlations would be reduced even further, 
with the EXVP correlations of the extreme velocities, the correlations ) 
4 for the Stokes velocities, Figs. 8, 9, are better for the crests than for the troughs i 
the waves considered. For the zero “up-crossing waves, the crest velocities 
for the elevation above the mean water level are biased on the artis and 
_ CREST VELOCITY VELOCITY , TROUGH VELOCITY 


@ 51 FEET 
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Bee, 
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NOTE: 1 FOOT PER SECOND = 0.305 METERS PER SECOND 


of OTS- Measured Extreme Horizontal Velocity with 


Velocity 3- 
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CALCULATED VELOCITY (FEETPERSECOND) 
FIG. 9. —Comparison of OTS-Measured Extreme Horizontal Velocity Magnitude with - 
tor Predicted Zero Down-Crossing Velocity (Waves 1, 3-10) S—— 
e for the | remaining elevations, on the high side; the net bias is minimal, but 
the scatter is about 20%. The trough velocities show - significant bias and 
considerable scatter—as much as 40% for the following trough. The same general 
Ss can be made for the zero down-crossing waves; however, for these 
waves (as compared to the zero up-crossing waves) the bias and scatter are 
"slightly larger for the crests, but smaller for the preceding troughs. ye enrol 


Although the Stokes-wave representation 


i 2° BIAS: 200% 

7 2 4 6 8 10. 

> wal 

“Stokes 

: 

| 

2 Ir BIAS: 20% 


earlier, that this ‘is only particular waves for which 
the | preceding trough | is larger than the following trough. For waves where the 
_ following trough is larger than the preceding trough, the up-crossing correlation | 


Therefore, for a large group of waves, the up- and down- 
a will not be distinctly different, and the scatter for each Stokes | 
‘representation w will be co comparable to the combined scatter shown here for both , 


a The EXVP procedure for calculating wave velocities can treat with accuracy 
4 single symmetric waves of relatively arbitrary crest-to-height ratios ; for a ra given 
wave height, period, and water depth, including Stokes-type waves of fixed 
_ This procedure may also be used to treat irregular waves. ‘Its accuracy for 
such applications has been checked against ocean wave measurements from B | 


_ the Ocean Test Structure experiment. V . Wave profiles and corresponding velocities 
_ predicted from EXVP for these waves show good agreement with the measure- 
_ ments for the entire wave period and water column. The EXVP crest velocities - 
were in closer agreement with the data than were the extreme trough velocities. _ 
_ The EXVP predictions matched the data significantly better than the common 
_ Stokes wave representations, especially under the wave troughs. -EXVP 
crest velocities showed negligible bias and less ‘than about (10% scatter; the 
extreme trough velocities exhibited no significant bias, but the scatter was 
- 20%-30%. The Stokes crest velocities were biased on the high side and showed | \ 
a 20% scatter, and the extreme trough velocities had considerable bias (also 
on the high side) and more than 50% scatter. 
Finally, this study has demonstrated that, in ‘general, improved accuracy in 
predictions of velocities for irregular ocean waves is obtained if the theoretical 
a7 formalism used accounts for a greater portion of the wave profile than the 
7 ere wave of interest. It has also shown that the EXVP procedure has the 
6 peony to account for more than one wave with acceptable accuracy. In 
contrast, common Stokes wave can at best match a 
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- matrix coefficient i in Fourier expansion of velocity potential; — 
| matrix coefficient in Fourier expansion of velocity potential; 


wil distance from bottom to mean water level; 
E, Fourier coefficient of velocity potential expansion in space; 


E. "error in dynamic boundary condition; = 


* = = error in kinematic boundary condition; 


EXVPD = Stokes’ wave reduction of EXVP; p. 


ee . Fourier coefficient of velocity potential expansion in space; 
acceleration of gravity; 
wave numbers for harmonic waves; y 
wavelengths for harmonic w waves; 


Gy 


least-squares sum of errors in n boundary conditions; 


time; 
wave periods of harmonic waves; 


horizontal component of water particle velocity; 
= vertical component of water particle ve velocity; 
horizontal distance; 


= vertical distance; eee 


wave surface elevation from mean water 


_ velocity potential function; an and 
= angular wave 


bscripts 
wave surface; 
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R ROUGH TURBULENT BOUNDARY LAYERS 


aa & The possibility of 3 giving a mathematically simple and fairly accurate description 
(of the oscillatory rough turbulent boundary layer deserves considerable interest, 
in view of the practical importance of this physical phenomenon. ‘Under natural — 
conditions, the bed boundary layer below sea waves is often rough turbulent, 

and reliable knowledge about the varying Siuid velocity near the bed and the 

associated bed shear stress is of primary importance when estimating, for instance, 
the magnitude of sediment transport. In the laboratory, a rough turbulent boundary © 

: = may conveniently be produced in water tunnel experiments above a rough 


~ Considerable work has so | so far been made to delineate | the essential features 
of this kind of boundary layer. On the experimental side, the oscillating water 
tunnel experiments of Jonsson (6) and Jonsson and Carlsen (8) play a prominent 
role, and the theory of the present paper will later be compared with the 
___ observations in their Test No. 1, which was the most accurate of their two i 


io series. Other experimental investigations, using different experimental methods, 
have been carried out by Kalkanis (10) and Horikawa and Watanabe (5). ar 
On the theoretical side, the most detailed mathematical treatment of the ‘rough 

_ boundary layer is probably the one of Kajiura (9). The tenor of his treatment | 

is to divide the whole boundary layer into three sublayers: an inner, an overlap, 

- and an outer layer, each layer with its characteristic mean turbulent viscosity. — 
It should be borne in mind that Kajiura’ s theory is a mean field theory in 
the sense that the turbulent viscosity v, = v,(z) is . taken as a function of the ‘ 
height z above the (theoretical) bed only; the time variation of v, has been 
ignored. Of course the omission of the time dependence of v, is an approximation 

_ for which the validity can -oaly be checked a posteriori by comparing with > 

accurate ¢ observations. Also, i it t should t be noted that Kajiura’s s ‘theory « deals = 


enough to o produce « an n overlap layer; f for saadiat excursion n amplitudes the overlap 


_ Now Kajiura’s theory with its separation of the whole boundary layer so | 
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= 
: three sublayers is formally compl complicated in that it require: requires ; matching of the analytic — 
solutions of the basic differential equation at the two boundaries between the 
_inner-overlap and the overlap-outer layers. Although the agreement between 
this theory and the Danish measurements indicates that the physical basis - 
; the theory is sound, one may wonder if the theory is not, after all, unduly * 
_ complicated. The present paper addresses itself to the following question: is 
it possible to reduce the mathematical complexity of Kajiura’s theory and yet — 
maintain reasonable good agreement with the observations? Specifically, the © 
: consequences of assuming a two-layer (instead of a three-layer) model will be © 


worked out. The change involves omitting the inner part of the boundary ad 
from the description. The advantage one obtains from this performance is that 
no matching procedure is to be carried out at the lower end of the —> 
{ layer. The price that has to be | paid for this simplification i is that the analytic a 
solution for the fluid velocity, in contrast to Kajiura’s solution, has to be 
terminated at a fixed level z = z, above the bottom. Such a situation is in _ 
fact well known from the theory of rough turbulent unidirectional flow with 7 
a a logarithmic velocity profile, the expression for the fluid velocity has no meaning 
. below the level at which it becomes equal to zero. . Moreover, ‘it seems to be 
rather commen among more recent workers in the field to choose a way of 
= approach where the analytic solution is terminated some distance above the 
_ The answer to the above question is affirmative. It is possible to reduce 
. _ the model to a two-layer model and still obtain reasonable results. The opportunity 
will be taken to discuss the position of the upper limit of the overlap layer, 
” Finally, the reader is referred to the recent. general treatment of Jonsson — ] 
on rough turbulent boundary layers (7). This extensive work focuses attention | 
on the existence of universal laws, i.e., the “‘law of the wall’’ and the ‘‘defect 
law for fluid velocity,’’ which seem to hold in boundary layers of this kind. 
Moreover, this work contains a survey over earlier developments in the field. A 
Assumrrions ano Equations Hee coral on On 


taal he present theory is assumed to hold in the region z => z,, where : 
0 defines ‘the theoretical bed level. If the bed is very rough, as it may be. 
when it is covered by artificial ripples, then the theoretical bed ievel will be 
lying somewhat lower than the top of the ripples. There seems at present to’ 
be no rule from which one can predict the magnitude of this small displacement 
theoretically. In practice, one determines the level z = 0 when fitting the observed 
velocity profile above the bed to a The quantity 2 = 
will be determined from the 


(l 


in which k , = Nikuradse’s parameter. k, can be determined experimentally =— 


by fitting the observed velocity profile to a logarithmic profile, as Shown 7 


Jonsson and Carlsen (8). If this distribution shall be useful, k, must obviously 

be a quantity independent of the wave phase. Probably, k, depends on the 

bed geometry only, so that it has the same value in the case of steady unidisoctional 


— 
| 
~ 
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ky: = in which + n= = . ripple height. gem new? | 
a = Consider now the region above z = . As mentioned above, there will in - 
be an overlap layer if the excursion mt 


form of the mean n turbulent viscosity v, in this layer. There si seems to be wither 
wide agreement among research workers ¢ 1a the linear distribution i 7 an a 


= z 


in which hich « = Karman’ constant. clear water, as assumed K 
= 0.40. Further, ui, = maximum of the bed f: viction velocity. Eq. 3 will henceforth — 
_be adopted in the overlap layer. 
_ The next point is to estimate the position of the boundaries of the overlap | od 
layer (9) estimated the lower to be located at the height 


= (Fig. 1) Making use of this observation, what will be dene here ie 
to assume the simplest imaginable picture, namely an overlap layer which extends — 
_ as far down as to z = z,. It seems natural to expect that the small uncertainties a 
in the location of the lower boundary should not significantly affect the magnitude 
of the nearbed fluid velocity or its phase. “de 
_ The overlapping layer will be assumed to extend to z = A , where A is as 
> yet undetermined. When comparing theory and experiment ame two different — 
"proposals for “a will be investigated. In the outer layer (z > A), it will be | 
assumed that v, is a constant. _ Moreover, continuity of vy, will be eines at 
=4, ‘so that de 


er. The constancy of v, in the outer layer and the continuity 
4 Hen are requirements similar to those made by Kajiura. At this point — 
& ‘one m might wonder i if it _ would not hi have been _more Physically reasonable to 
assume a form of v, in the outer layer : that decreases with height. For instance, — 
= ‘ Lundgren (11) assumed a form for v, that increased linearly with z for small | 
; 3 z and decreased exponentially with z for great z. A form of this kind may 
~ actually be more physically satisfactory; in consequence, the basic equation 
of motion, Eq. 9, in general becomes solvable only by numerical methods. ; 
a In the outer layer, there is great scatter in the experimental data. Probably, © 
the idea of introducing a time independent turbulent viscosity is not very adequate 


here. It seems, therefore, to be to keep the form of v, as 


h 
| 
a —— satisfied. In the following, it will be assumed that the physical conditions - 
{ 
| 
- 


di 


above the boundary layer, +(z,t) = shear stress, and p= is of water. 


cuita 


sonia: 


Value of in cm? 
FIG. 1 Viscosity. ‘Data from Jonsson and Carlsen (8) 
In Eq. 5, it has been assumed that the nonlinear ccnvective term is negligible. 
it has been assumed that the boundary is 80 thin that hat the 


from potential theory. 


turn now tothe basic § 


neasuring the deviation the free stream velocity: 
The idea of a - position dependent but time independent hol means that 


It is now to introduce complex notation. stream velocity 
a be written as U(t) = U exp (iwt), in which U = real velocity ages 


= which u(z) and u,(z) are in . general complex. due to the phase shift relative 
free stream velocity. Eq. 5 may now be written 


It is fortunate that this oul is ‘explicitly ‘solvable both when v,(z) varies ag 
linearly with z and when it is constant. The given in the 

ve 

= Consider first the overlap layer i in which v on is _ by Eq. 3. Invading 
the nondimensional variable, 


the nondime | 


urns out “that Eq. 9 reduces to the differential for the 


in which D is a constant quantity defined by Ms 


+ ker €,) (ker, — keié,) beié,] 
In this expression, bei etc. = first order Kelvin functions (1). The 
expressions for &, and €, are given by Eq. 10 with z = z, and z= A, respectively. 
= 4 ‘The quantities A and B are nondimensional constants, in general complex, that y 
er have to be determined from the boundary conditions. The reason it is ae 


| 
— Onvenie O write the solution ~ 
4 
- 


If, by ‘contrast, the ¢ linear at v, z were 
assumed to hold for great values of z, then one might put A = 0 because sid 
g the magnitudes of ber and bei increase rapidly for great values of é. ae 
Now consider the solution of Eq. 9 in the outer layer, in which v, is given 


in which .....- . 
Further, T is a new constant, in general complex. In writing Eq. 13, use has 
‘been made of the requirement that u, sade 0 when z 23m, The minus sign in 
front of the equation has been added for convenience, as the physical defect 
velocity is ordinarily negative. The solution for u, in the outer layer may of 

course also be expressed in terms of » but it is more simple to express it 


directly i in terms of z, in Eq. 13. wan estan o 
It is, in principle, straightforward to determine the complex constants rm 
and T from the bey conditions at the two interfaces of the boundary 


A A (ber 4 ibei€, B(keré, + ikei € , 
‘Second, continuity © of uy at =A implies t that rom 


Third, there is the condition of continuity of stress at = As 
v, has been assumed continuous at this level, it follows from Eqs. 7 and 8 


r: that du 4(2)/dz itself is also continuous. This condition leads to the equation © > om 
A (ber, + wre ‘Biker, + + ikei = iT 


Here use = been made of the relations (1) 


Eqs. 15, 16, and 17 the constants may be determined: 


A = ker,&, - 


kei + 
ber bei €, i(be ei 7 eré,) a 


ert, kei + bei,€, ker €, ker,€, beig, — beré, 
The relative simplicity of the pane for A and B in Eqs. 19 ‘and 20 was 


it is, in 1 principle, Seapets es to calculate D, A, B. and T from Eqs. 12, 
and 19-21 and thereafter calculate u,. The Kelvin functions may in general 


| 
| 
| 
| 
oa 
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be found from tabular works (1). ‘that in the cas case 0 se of small 
for instance, when & = &,,, the approximate analytic expressions to the Kelvin — 
functions are often useful. These expressions will not be written down here. 
‘It may be worthwhile to give the formal ee for the magnitude “E 
phase of the defect velocity: Oars 


in which i‘ otdiventions follow Jonsson (7) (therewith , = 0 at the point 
in which the fluid velocity u vanishes). Now define, in Fae eet notation, 


=-- 
M cos (o, — 6, + ©) + cos 
& 


; In the expressions on the right, the only quantities varying with z are M, _ 


8 and 6. and of the defect velocity are now calculable 


= constant; “the linearity 
~such has no direct bearing on the fact that the. present boundary layer is a 7. 
turbulent one. Qualitatively, the aan al variation would be found in a laminar 


for all values of z. 


. 
_ Now compare the above theory with Test No. 1 of the Danish.2 measurements | 
6, 7 - In this case, the existence of the overlap layer was obvious, as A,/ky __ 


q 
The final expressions then take the foll a 
— 
7 
7 
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see Eq. 2. The wave | period was T = 8. 395, the 
ky = 0.91 in. (2.3 cm), and the free stream velocity U = 83.1 in./sec (211 — 
_em/s). The theoretical bed level was Gaplaced 0.098 in. (0.25 cm) downwards 
relative to the top of the bed ripples. When fitting the experimental velocity _ 
4 profile for u(z,t) to a logarithmic form, Jonsson and Carisen (8) obtained values , 


for p at at the bed ranging from 62.0 in. /sec* (400 cm?/s? )-—72.1 in. /sec? - 
@ —465 cm dim ). This indicates the following choice of friction velocity: u, = a 


,/p)'? = (430/p)" '/? om/s = 20.7 cm/s (8 (8. 15i in. /sec), which will henceforth 
be adopted. Correspondingly, 0.167, 4 
Using the experimentally determined values for u(z,¢) and t(z,¢), one may 
calculate” time dependent values for v v,(z,¢). For. convenience, 
values from Jonsson and Carlsen (8) have been plotted in Fig. 1. The points - 
oh to measurements differing 45° in phase. Negative values of v, have been | s 
_ omitted. The points in the figure indicate that it really makes 4 Metals 
_ approximation to represent v, by a time independent, linearly varying function 4 : 
Be of zin the overlap layer, as was is anticipated a above i in writing Eq. 7. The straight 
is in the lower part of the figure represents Eq. 3. The inner and outer limits 
the overlap layer are, as one must not fixed very 


| 


— 


FIG. 2. Defect Relative to to i, = 20.7 cm/s, 8 = Gem. Data 
from Jonsson (7), Test No. 1. —: A= 1/28; —--:A= (Im=3.28ff) 
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the data. _— seems to te no experimental evidence against extending the 
overlap layer down to z = z,. Greater uncertainty is associated with determining 


- the outer limit of the overlap layer. To investigate t to some extent t the ‘sensitivity 


ti. 


4 


FIG. 3.—Phase of Detect Velocity — 6 cm). Data from Jonsson (7), Test No. 1. a 


for the magnitude of A will be worked out here. The first one isto put 
1 Tig, thas choice to give af 


(Corresponding t to the vertical line to the right in Fig. 1), in which 8 = boundary 


| 
| 
&g 


‘layer t as defined by Jonsson (7,8). In Test No. 1, was determined 7 
to be equal to 2.4 in. (6 cm). This case corresponds to&, = 1.042. havo ‘a 
As the second for 4, Kajiura’ s (9) will be invoked: 


in which @=0.05 | 
in whic i 


e = friction factor. In Test 1, one calculates d= 0. 33 i in. . (14. 1 mm). Correspond- - 
ingly, = 714. This case is also i in | l. a 


ace 


- 
ionless Fluid Velocity, Relative to ° 4, = 20.7 cm/s 6 = = 6 cm. ‘Data 
| 8), ‘Test | No. 1. 
‘4 Figs. 2 and 3 show calculated amplitudes and phases of the defect a 
u,, and Figs. 4 and 5 show corresponding quantities for the fluid velocity uw. 
The from Test are also seems in fact 


fonction 


- that the two- layer model investigated here yields results i in agreement with the 
observations. The gross features of the empirical variation of amplitudes and om 

7 phases are reproduced by both the two input values for A. It is somewhat 
surprising to see how insensitive the final curves are with respect to A; + 
Fig. 1 the difference between the turbulent viscosities in the outer ‘region = 
the two cases is actually quite great. Which of the two alternatives for A yields = 
the best niennyyy with observation is not answered definitely by the ape 


10 


FIG. 5.—Phase of Fluid Velocity (5 = 6 Jonsson (7), Test No. 1. 


ag in the region close to the wall (z < 8), it seems yet as if the simple choice, 
be A = 1/2 3, fits the data most accurately. Especially regarding the phase of 
; the fluid velocity (Fig. 5), this choice for A appears to give best agreement. bom 
- Kajiura’ s (9) gives his theoretical amplitudes and phases for the fluid velocity 
; wt under the physical conditions | of Test No. 1. His figures are thus directly _ 
ve comparable to Figs. 4 and 5 in this paper. Inspection of the figures shows — 


| | 
— 
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that the agreement between ‘theory is almost same as that 
in the present case (with A = 1/2 8) as in Kajiura’s case. ging MAL cet 

a. Finally it ought to be mentioned that the -urves in Fig. 3 are not very different — 
from the straight line determined by the following formula, proposed by Jonsson — 


This formula \ was found to agree with several experimental tests, it, 
warrants considerable interest from a practical viewpoint. 


1. The main purpose of the paper has been to develop a two-layer theory | : 
the oscillatory rough ‘turbulent layer, thus, aiming at a reduction of the 
mathematical complexity of Kajiura’s three- -layer theory (9). In the general c case 
of large excursion the total layer was assumed to consist 


/ from z z=A embents. Two different choices wore made for the > input parameter — 


2. The calculati on rests. on that the periodic variation of the 
exp (iwt). Specifically, this assumption was used in deriving Eq. 9 from the 
basic ay of motion, Eq. 5. Experimental evidence, such as the Danish 


function in | time. Thus, rem the fluid — sul the shear stress are in n the _ 


‘is, after all, basically nonlinear, and this particular behavior encountered here 
is a consequence of the simplifying assumptions being made: first, that the 
nonlinear convective is omitted in basic ‘Eq. 5; ‘second, that t the time 


in 1 this ¢ context to refer also to the discussion by Jonsson (7) on the influence 
the nonlinearity of the problem has on the proposed “‘laws of the wall. oie 
_ 3. The results of the present theory are given for ‘the two sublayers in Egs. 
_ 10-30, and the results are shown graphically in Figs. 2-5 for parameter values — 
_ pertinent to the Danish Test No. | (6,7,8). Reasonable agreement with 4 
observations is in fact found, thus « confirming the usefulness of the two- layer 
model. In the nearbed region, it appears that the choice A = 1/2 8 (Eq. 31) 
_ gives accurate results. It should be noted, though, that the results are not very 
4 _ critical with respect to different input values for A. In conclusion, the obtained 
results. are re satisfactory, not ot only because they show the usefulness of the ‘two- layer 


the mathematical treatment of the problem. ots ath 


| 
4 
> q 
— 
| 
they generally support the idea of introducing a turbulent viscosity, which is 
n. This feature facilitates 
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followin sy mbols are used in this pape 
quantity given in Eq. 19, i 
ball nearbed excursion aaa of fluid particles; i. 
= quantity given in Eq. 20, B = amplitude; 
= quantity given in Eg. 12, D= amplitude; en, 
= Kajiura’s outer limit, Eq. 33; pray as 


= friction factor; 


= Nikuradse parameter; 


quantity given in Eq. 21, T= amplitude; 


velocity bo free s stream; 
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1 
ere velocity, sity, a, = = ‘amplitude; 
= height above theoretical bed level; dation ant 
height of outer limit of overlap layer; 
phase defined in Eq. 23; 
phase defined in Eq. 23; augas? 
quantity defined in Eq. 10, €, and é referring to z = z, and z 


= phase of u,, Eq. 22; 
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= By Todd ‘Walton, Jr., M. “ASCE Francis F. ASCE 
4) 
_ Various investigators (4), Van de Kreeke (8), King (5), Goatide (2), Ozsoy 
(), Es Escoffier and Walton m (1)] have proposed solutions to the hydraulic flow ~ 
a - through inlets problem m for : simplified type inlet bay systems in which the ocean 
: _ tide has been assumed sinusoidal and both inlet cross section and bay surface ; 
area are constant. Of these solutions, only Ozsoy (5), and Escoffier and Walton a 
©, have retained the inertial term and provided an analytical solution . Addition- — 
ally, Escoffier and Walton (6) considered the case of "tributary inflow in their - 
"analytical approach. Ozsoy (5) solves the ‘inlet- bay | problem with inertia using 
“the describing function technique frequently used in the theory of nonlinear 
control systems (7), while Escoffier and Walton (6) use a “least squares” , 
linearization approach tothe problem. q 


‘The present note describes a solution technique using a linearization technique = | 


of equating friction- work over the tidal cycle in a ‘nonlinear approach to 


_ friction-work in a linear approach, and thus avoids apriori assumptions of the 
- form of bay tide. This method was first used by Lorentz (8) in tidal computations _ 
7 of the Dutch delta works and more recently has been applied by Van de Kreeke © : 


n of motion for one dimensional channel flow 


which the advection of momentum term hes been and u = inlet 
_ Channel velocity; g = acceleration of gravity; h = water surface elevation — 
channel bed; R = channel hydraulic radius; and f = Darcy-Weisbach friction 
- factor. The second term on | the right hand side of Eq. 1 can e linearized — 
the method of Lorentz (8) in the following equation 


a 
tidal work non-linear = tidal linear. 
Engr., Evaluation Branch, Coastal Engrg. Research Center, Kingman an Bidg., 
Retired; formerly Chf., Hydr. and Sanitary Section, United States Army Corps of 
7 Note. —Discussion open until January 1, 1982. To extend the closing date one month a 
a written request must be filed with the Manager of Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication on March 26, 1980. 
_ This paper is part of the Journal of the Waterway, Port, Coastal and Ocean Division, — 
: Proceedings of the American Society of Civil Engineers, ©ASCE, Vol. 107, , No. WW3, 


LINEARIZED SOLUTION TO INLET EQUATION 
= 
— 


in which U,, = maximum channel velocity. 


Integrating from from the ocean ead of of the to the 


u+h 


in which F = f{L/4R = friction parameter; h, = elevation of 


Combining Eqs. 4 and 5, assuming a ocean using the 


a 


in which o = T = tidal period; and a, ocean tide 
A dimensionless linear differential equation ‘can be 


solution to equation 6 can be to be RAE 

in which {= Ja, = factor or ‘‘amplification”’ “factor; 


water surface in bay; 4, = elevation of water surface i in ocean; and L= effective 4 


in which cross- area in gorge of inlet, = water surface area 


= 


monic channel velocity and performs the requir - 

) 
 &§ 

> a 

| 

— | 
fae 
: 


7 ‘As B has an implicit measure bay response in it, equation 
7 to find U,, and thus ¢ and ¢ in explicit form. An equation for U,, can be 

_ found by equating the tidal prism over a half tidal cycle to the flow through 


inlet in the same time period. we 


Tidal prism = 2a, A = (10) 


Eq. 10 U,, _can be found as 


Sy ‘ata 8, 9, au 11, the following explicit representations of bay tide See 


{, and bay lag € can be found 


¢ 
4, 
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? of Escoffier-Walton (1). 
Graphical solutions of Eqs. and 13 in ‘Figs. | and 2 for 
‘various values of the ‘ ‘damping”’ coefficient. It is noted that all solutions presented 
_ in Figs. 1 and 2 correspond well, with three of the solutions providing exactly y a 
“the same answers ©), ‘Walton on. and present solution]. 
fg ‘Orsoy (1077) 
King (1974) | 
g 


q 


= Dimensionless Tidal Frequency 


= 
» 


= Keulegan ‘ ‘repletion”’ coefficient . 
"a preceding equations and Figs. 1 and 2 are not plotted in Escoffier-Walton - 


O) but are detailed here to add further me insight into the problem and 
justification for the method of 


ACKNOWLEDGMENTS 


want conducted at the Coastal Engineering Research Center under the Coastal 


; = Research Program of the United States Army Corps of Engineers. _ 
Permission to publish this information is 


be 
| 
— 
‘ 
nd 23 of that — 
4 
| 
a 
4 


F. F., and Walton, T. Li, ,Jr., “Inlet Stability Solutions for Tributary Inflow,” 
be Journal of the Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 105, No. 
WW4, Proc. Paper 14964, Nov., 1979, pp. 341-355. aw 
2. Goodwin, C. R., “Estuarine Tidal Hydraulics,’ thesis presented to Oregon State 
(7 University, at Corvallis, Oreg., in 1974, in partial fulfillment of the requirements for = 


3. Graham, D., and McRuer, D., Analysis of Nonlinear Control Systems, Dover Publica- 


; 4. Keulegan, G. H., ‘Tidal Flow in Entrances Water-Level Fluctuations of Basins in 
a Communications with Seas,”’ Technical Bulletin No. 14, Committee on Tidal Hydraulics, Rad 


Army Bagineer Waterways Ex Station, » July, 


5. King, D. B., “The Dynamics of Inlets and Bays,” Tetincleg) Report No. 22, University — 
Lorentz, H. A., *‘Verslag Staatscommissie Z Zuiderzee 1918-1926, Staatsdrukkerij, The The 


c 7. Le Van de Kreeke, J., “Water Level Fluctuations and Flow in Tidal Inlets,’’ Journal 


4 7 of the Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 93, No. WW4, Proc. | 
108 Paper 5575, Nov., 1967, pp. 97-106. 


ILITY or R RUBBLE Mounb BREAKWATERS 


Tedd Wailea, Jr.,' and J. Richard Weggel,” Members, “ASCE 
Numerous investigators [Ahrens and McCartney (2), Thomsen, et al. (1972), 
Bruun, et al. (2), and Losada and Gimenez- Curto (6)] have shown that wave — 
steepness, and consequently \ wave period, is a factor in the stability of rubble — 
mound breakwaters. Some present design methods (5) do not explicitly consider — 
wave period. A modified wave steepness dependent stability number, accounting» a 
q for inertia effects is presented herein. The following derivation of an equation a 
for armor stone weight at incipient instability (motion) generally parallels an 
earlier derivation by Raichlen (6), but retains lift and inertia terms in the final 
analysis. As a result, a stability equation similar to that of Iribarren and Nogales 
gs is found, but, which includes an inertial parameter dependent on wave 


'Hydr. Engr., Evaluation Branch, Coastal Engrg. Research Center, Kingman Bldg., 

"Unites States Army, Ft. Belvoir, 
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Rubble mound —_ stability is reduced when the inertia 
_ Raichlen (6) has shown that when moments are taken about the point < ol 
‘contact (support) of an armor unit for the condition of incipient motion as 
‘shown in Fig. 1, the following equation canbe found: = a 
(Fp + F,)c,lcos® + F,c,/sin® = —a) 


& in which =— Coles P 


=. (c; P 


‘and ¢ and = coefficient, the lift and the 


__-virtual- mass coefficient ‘respectively; @ = the angle of repose of the armor 


Point of Rubb 


ion of Armor Units [from F Raichlen (6)} 


unit; the angle the structur Slope makes with the horizontal; 


characteristic length of the armor unit; and c i= = the distance from the point 
of support to the center of gravity of the unit; cl’ = the ‘projected cross 
_ section area of the unit on a plane perpendicular to the direction of the velocity; 4 


= the volume of the unit. 


We Ps ees 


own assuming that velocity * me slope is positive that 


| 
| 
) | 
4 
4 
1.—Definition Sketch for Incipien 
4 
. 4 &g 
q 
and tha ‘ 
| 
and 
| can b | 


a shallow water theory the velocity can be written Y 


= water surface elevation above the mean water level; C = one 
eelerity = V gh; andh=waterdepth inca 


The water surface elevation is by Wis A 


eh A 


0) 


{i 2 A end 
in w which height; = wave angle = ot; o = 2n, 


phase time; and T= wave 


(Gye - oo A — 2) 
— 1)(tan®@ cosa — sina) + Cy eo cos ) BA 
J 


‘Assuming 1 no lift force (C, = 0) and no inertia force | ej = ; 0) an | equation 
similar to Iribarren and Nogales’ ay can nSs found if the peak aoe velocity 4 4 


4)? (tan 8 cos a — sin a) as WE In 

By ‘retaining inertia and lift terms, and ‘shifting the ‘wave phase angle, 


—ulu| 1+ —— 
. &§ alt du 
pie. 
 inwhich .. 
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W, by to reflect ‘downslope velocities an and Eq. 16 can 


(145 ) sin? [sin 
sin” sin’ 


- 8 cos sin a) — Cut (7 
a 7 pr that the inertia force can lead to an instability phenomenon “S 
which wave steepness is a critical factor. 
- As present model tests of rubble mound structures do not explicitly consider 
- _ inertia, , stability numbers of the ‘‘Hudson” (9) type such as N or K, (in 4 
} Ky, = (tan 6 cos a — sin «a)/N cot a should be looked at critically when ~ 
; rubble mound siructures are to be built utilizing laboratory stability numbers 
based on experiments conducted with lower than those actually 
On dividing Eq. | by (S (tan @ cos a a modified 


q 


tan 


-inwhich W’=- RET 
(S— 1)’ (tan cos — sina)’ 


C, 


can be thought of asa‘ ‘phase stability” number and R, 


As typical design would be in accordance with an equation of the type (oak aie, 
wea give a rock weight equal to W’, the designer should insure that site| 
a conditions do not produce a maximum P, value than that 
in the laboratory tests on which the stability is based. 
As an example, typical rubble mound breakwaters encompass ‘slope angles 

from 1 on 3 (18.4°)-1 on 1.5 (33. 7°). Typical design wave steepness might 
ange from 0.01-0.10 (in which 0.14 is the deep water limiting wave steepness). 
Assuming an angle of repose, tan 0 =_ ‘1.0 [see, e.g., Hedar (10)]; S = 2.65; 
and an inertia coefficient Cy = 1.5 (for a hemisphere on a plane surface); 
al typical range cof R values: would be R = 0.045-1. 030, with the higher va value 
to steoper slope (1 on 1.5) and larger wave steepness (H/L 


Plots of W max angle corresponding to m maximum versus R, and 
7 P, versus fad are given in Figs. 2 and 3. These figures provide some “insight a 
ns for which inertia effects become important. ; 


| 
i \ 
| 
4 
4 
\ 
| 


thet: as R increases increases until as R 1.0, P, is: much 
larger than it is for the case when inertia conditions are not considered (R- 

= 0, P, = 1.0). It is also apparent that high wave steepness and large subtlomcund 

slope angles could be critical in design should laboratory model tests not 
encompass the same range of wave steepnesses characteristic of the site. 

;  Egs. 20 and 21 may also be used to establish a maximum allowable structure 


slope . Structures bulls having slopes steeper than the allowable slope run aa 


iv ot 


ae FIG. 2 —Plot ‘of Phase Maximum Inertia- Parameter Ratio 


“FIG. 3. —Plot o of P Phase Stability Number Inertia- -Gravity Parameter Ratio 


risk of damage due to inertial effects. Fig. 3 suggests that values s of R greater : 
- about 0.5 will 1 result i in values of P, greater than about l. 25 ¢ or r the required 


‘maximum value of R can be established beyond which an unacceptable risk 
of failure exists. The maximum wave steepness H/L and inertia coefficient, : 
Cy, and the value of S along with the angle of repose for the | armor, 6, then 


een is used. For larger values of R, the situation is worse; +; thus, “some q 


=f 


| 

| 
| 
| 
x 
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‘determine a maximum structure slope. From Eq. 2 
i 
‘tan 60s a — sin a = 
in which the value of H/L is the largest teiia waves steepness at _ 
site and R is the largest acceptable value based on the largest acceptable P.. . 
Eq. 22 is presented graphically on Fig. 4. ‘The value of a equation 


30 
Se FIG. 4. for Steep st 
22 slowable slope. For ex example, if 
(H/L) = 0.05, S = 2.65, tan 6 = and the maximum R is 0.5 
(corresponding to P, = 1. 25) then 


Fig. for 0, = 33 33. @=1. 51. 


4 

> 


7 ks 1.5 could thus experience pia if the stability tests upon which the design 
based did not | values H/La as great as H/L= = 0. 
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mean 1 velocity over a wave at any point’ 
- troughs is exactly zero, thus causing zero mean mass flux across any vertical 
control section between the bottom and a fixed point deeper than the troughs. — 
a Second, the Lagrangian mean velocity of all water particles is a nonzero second — 7 
7 order quantity, and is always in the wave propagation direction. At first sight, 
it seems justified to ask, as probably many students do, how can it be that 7 


- control section remains zero? In the present note we ere , what is hoped — 


to be, a clear and descriptive answer to this question. 
i In Fig. As we 5 show the pathlines of four selected water particles Going one 
wa 
section AB: are » moving upwards. Particle 1 crosses the section AB twice during 
_ @ wave period and contributes nothing to the mean mass flux. Particle 2, which al 
starts from a location nearer to AB crosses the control section only once, from 
left to right, thus | contributing to a positive mass flux. (The writers define the 
in the wave propagation direction positive. All particles in zone 
1, (the zone shaded with horizontal lines) are the same type : as , particle 2. The 
width c of this zone is identical with the particle displacement during a wave 
period and is given, according to Phillips, as follows: 


‘ which a is the wave ‘eeteb, k the wave-number and A the mean water 


depth, 


The area of zone 1, , is given by 


wn ‘can be shown , accurate to order a’ for the mass flux, that all particles: 


a “Research Engr., Coastal and Marine Research Inst., Technion-City, Haifa, Israel. 
Senior Lect., Faculty of Civ. Engrg., Technion-Israel Inst. of Tech., Haifa, Israel. = : 
_ Note.—Discussion open until January 1, 1982. To extend the closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, 
_ ASCE. Man ‘uscript was submitted for review for possible publication on March 11, , 1980. 
This paper is part of the Journal of the Waterway, Port, Coastal and Ocean Division, 
Proceedings of the American Society of Civil Engineers, ©ASCE, Vol. ny No. WW3, 
1981. ISSN 0148- -9895 /81 /0003-0202 /$01. 00. 


| 
— 
while all individual natticles move forward the mean mase flux across the 
| 
J 
where H is 
: 
| / 
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; ‘fedide zone 2 (the ellipse shaded with vertical lines), which has semi- axes a 
‘ and b, have pathlines similar to that of particle 3 and cross AB only once, - 
right | to left. All these particles contribute to a negative mass flux. i 


as by Phil Phillips. 2m od} 4 leubrvibat its 


ae we obt in the are area of zone 2 (hich is is denoted asA,): — com 


2 sh? (kh) 


: Particles, in the neighborhood of zone 2, that start above z = —H (e.g. particle 
4 a , do not cross AB at all. Particles, in the neighborhood of of the he cllipes, that 
start below z = --H behave like particle t. 


FIG. 1.—Pathlines of Typical Water Particles Waves 


9 


From Egs. - 2 a and 4, we see that the mean mass flux during a wave period, 
given by (A, - A,)/T in which T is the wave period, equals zero, despite 
fact that all particles move to the right. 
. hh order to avoid misunderstandings, note that the choice of any | different 7 
initial instant, would result in the same discussion and results. The only difference - 
= will be ‘the displacement | to ‘the right, bya distance not bigger than 2a, of — 
the ellipse (zone 2, in Fig. 1). It is worth 1 mentioning ‘that, when point B is 
located in the region z(—a,a), the elliptic zone 2, is partially occupied 
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at the initial instant. Thus the area of the relevant part of zone 2 is decreasing © om 


_ from ab, for B at z < a, to zero, for B at z = a, causing the mean flex 
a°/T coth(kh ), as should be expected. tight mov! 


to increase from zero to 7 
j A new approach is ‘ied in this short note to show that, despite the fact 4 
that all individual water particles move forward, the mean mass flux — 


a vertical control section (from the saieech to a poms located er wave 
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Mixine or Buoyant Surrace Jet over SLopinc Borrom* 


Closure by Bijan Safaie,” A. M Ye 
would like to thank Adams for his interest in the writer’s work 
ad for his valuable discussion. It appears further clarification and discussion — 
> on the mechanism of separation and the corresponding normalization is warranted. — 
; This, together with the writer’s comments on the points brought up by the © 
discusser, is presented i in the following. was Lot om tee 
_ The writer’s observation, from the side glass window, indicated that the — 
‘mechanism of separation was very similar to that of a buoyant ball entering - 
a fluid with certain initial kinetic energy. In fact, one can show that, neglecting + 
friction, such a ball can eee in the vertical direction to a limited sonae, 
say L a V/s’, in which L = the maximum penetration . of the ball; V = © 
| the velocity of the ball in the vented direction; and g’ = the modified gravitational 
acceleration based on the density of the ball and ‘the fluid. The movement 
of large-scale turbulent eddies were very similar to that of the buoyant ball 
_ previously mentioned. This observation led the writer to select kinetic energy 
flux and ‘buoyancy flux as important parameters in forming the proper length — 


: ‘scale, L,. Based on these observations, the water depth at the separation point 
can physically be interpreted as the maximum depth which large-scale eddies 


can penetrate. Based on these comments and one available ee data, 


x The normalized maximum penetration om isa function of source den- . 
etric > Froude number, F.(h), only and is independent of both: aspect ratio 
- bottom slope. This can be seen in Fig. 11 where data points correspond — 
to various values | of aspect ratio and bottom slope. The same functional — 
relationship is expected for maximum penetration depth of surface buoyant 
h: jets discharged into > deep water. Carter’s data shown in Fig. 11 substantiates | 
‘this hypothesis by noting that his data is for BSJ discharged into deep water 
and that they follow other data points closely. Therefore, experimental evidence © 
_ does not support the finding in Ref. 24 that = constant in Eq. 33 for sare 4 
> Although the dimensional analysis and normalization of h, presented =i 
the discusser is valid, the normalization which was used to obtain Fig. 11 has" 
an advantage in unifying the data points which were obtained under dissimilar 
ae In general, based on the experience with submerged buoyant jets, Q, and 
A, were expected to have small effects on the dynamics of the BSJ, however, e 
is not a valid assumption for BSJ over sloping bottom. Note that 
“November, 1979, by Bijan Safaie (Proc. Paper 14954), 
a Asst. Prof. of Civ.  Engreg., State Univ. of New York, Par ker Engre. Building, Buffalo, 
ONLY. 14214. 
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evidence ‘from various sources this point (see Fig. 13). 

_ The writer suggested that the lateral extent of the attachment region ~ 

: f be approximately equivalent to the nonbuoyant spread of a turbulent jet with © 

_ the same source geometry and velocity. This was based on the flow visualization — 
\g . using dye crystals and photographs. Other methods of determination might result > 
_ in different values, moreover, very accurate determination of the lateral extent a 


. of attachment region experimentally is rather difficult by noting that the angle — 


of spread for nonbuoyant jets is approx 9°-12°, 


Inter SOLUTIONS FoR TriBuTARY INFLOW" 
wil Closure by Francis F. Escoffier,” F. ASCE 


and and Walton, , Jr. 


_ The writers would like to thank Bruun and Mehta for their interesting discussions - 
; end for drawing attention to the notions of inlet stability as opposed to inlet 
= Bruun discusses the subject of inlet equilibrium where sedimentary 


“equilibrium” implies that an inlet channel - flushes as much sand away as it 
takes i in by littoral drift feeder currents. The writers refer to hydraulic ‘ ‘stability”’ 
rather than sedimentary “‘equilibrium”’ as the title of the paper implies. The 
writers concur with Bruun that littoral drift imposes different constraints in 
- the inlet. This has been taken into consideration by a varying factor p as given 
in Eq. 46 which depends on the method in which the channel f fills. It is recognized — 
that due to varying wave action on 1 different shores that different values of 
would resui. ‘edimentary equilibrium can best be expressed as Bruun (11) 
has noted by a v...stant inlet channel velocity for a particular exposed shoreline. _ 
_ It is mistakenly noted by Bruun that Walton and Adams (16) work on outer 
inlet shoals versus tidal prism does not account for wave action when, in fact, 
a criteria for wave action has been taken into account | by a wave energy ove ; 
(in which H = wave height; and T= = wave period). 
‘The writers concur with Bruun in his comment that ‘“‘it is dangerous to ae a 
fixed ratios or too strict generalizations on inlet behavior.”” As is too often 
the case, inlets seem to have a ‘“‘sixth sense’ of their own which continues | 
to defy the engineers attempts to quantify them. If progress is. too be made 
though, such attempts must be tried | and judged accordingly as ‘qualitative 
apne lack the substance on which | to base sound engineering decisions. mnt 
a “November, 1979, by Francis F. Escoffier and Todd L. Wakon, Jr. (Proc. Paper 14964). — 
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DISCUSSION 


Force Anatysis: AN ALTERNATE ProcepurRE* Vel. 


" ‘The writer would like to thank Garrison for his constructive comments. The 
writer fully realizes the complexity of the fl flow phenomenon and the limitation — 
of the method suggested in ‘the | original paper, and» wishes to clarify a few 
points that might have mislead the readers. 3 a 
_ The proposed method of wave force analysis was in fact based on the premise — 
that the forces acting on a section of a cylindrical pile can be described by © 
7 _ the particle velocity, acceleration, pile diameter and fluid properties (Eq. 1). 
4 The Physical meaning of this assumption was adequately described in the original — 
paper. The writer fully agrees with Garrison that the particle motion history 
_ is an important consideration. The proposed method did not neglect the effect — 
of motion history entirely, since the motion history is somewhat represented — 
by u, a and J, where / is directly related to the ‘‘displacement ratio”’ cited 
® Garrison. It is fair to say that a more complete representation of of motion — 
a history is needed for further advancing the state-of-the-art on wave force 
: calculation. However, the lack of a detailed representation on particle motion 


history (which is common in most research works on the related subject) should 
not overshadow the key points taken in the paper. It was not the intention 
a this paper to correlate the drag and inertia coefficients with the Reynolds | 
number and Iversen’s Ss Modulus based on instantaneous fluid kinematics. The 
“suggestion of the paper was to bypass the step of using the Morison’s equation 
2 altogether. Eqs. 7 and 8 were used to illustrate the fact that the proposed | 
' method is more ‘general than the Morison’s approach. In other words, if one 
accepts the use of Morison’s equation with a given pair of C, and C,, Eq. 4 
7 should indicate that the dimensionless force $,(t) is a linear function: of J 
with C,, being the slope and C, the intercepts. The relationship between 6, 
and J, as shown in Eq. 8, constitutes a hyperbolic function. This suggests that 
a plot of , or , versus J for various R can be made to calculate the C, 
and C,, if Morison’s equation is applicable. The intent of this discussion is 
to demonstrate that the Morison’s equation f. alls within the domain of the proposed 
‘method. But, the purpose of the paper was to encourage the use of the alternate — 
—— as proposed, which is to analyze wave force data based on Eqs. 
5 and 6 without using C,, C,,, and Morison’s equation. The proposed approach 
simplifies the data analysis procedure because the direct measurement of wave 
_ kinematics may be used to correlate with wave force data without relying upon 
any wave theories. In addition, the results of simulation may be applied to’ 
- design problems in a more ‘direct manner. | . For example, in many instances the 
_ water particle kinematics may be measured rather easily during : an offshore 


February, 1980, by Yuan Jen (Proc. Paper 15166), 
*Dir. of Ocean Engrg., Dr. Maxwell C. Cheung & Assoc. , 2081 Business Center Drive, 

Suite 110, Irvine, Calif.92715. Yo . inl! 


q 


operation while the measurement of waves is more difficult. The 


_ alternate procedure would allow wave force calculations to be made directly 
- from the measured wave kinematic data and subsequently applied to the design a 


- _ The existence of a formal second-order solution may not be in contention, 
as Isaacson states, but his original article (1) appeared to imply that it — 
of doubtful use because of the line of singularities. Ir In fact, the crucial point — 
does not concern just second- order solutions but the whole procedure, s starting 
with first-order solutions, for each order yields a linear problem like the one © 
discussed in the original paper. Only under special circumstances can one avoid 
these logarithmic singularities, and higher-order approximations cannot cancel 
those that occur "earlier. ‘Situations where the latter seems to occur are mot 


. that finite values for force and moment result. Similarly, the leading-edge 
Singularity in thin-wing theory does not prevent it from giving useful results a 
for the lift even though the pressure distribution is incorrect near the leading 

~The p possible occurrence in the linearized problem of a logarithmic: ‘Singularity — 

at the intersection of the equilibrium free surface and boundary walls was already — 
pointed out in 1910 by J. Hadamard (8), and further elaborated by G. Bovligand Pa 
(9, 10). In his second paper, Demgnan claims to prove that there are aso 

‘Singularities in the “‘exact’ ’ problem, but the writer finds the proof obscure. 

In the exact problem for a flapper v wave-maker, there must, however, be a 

: singularity at the intersection of the flapper « and the bottom, just as in the 

a Since submission of the writer’s note the logarithmic character of the singularity 
has been confirmed by Miloh (11) and a solution to the second-order problem 
Hunt and R. E. (12, 13). 
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10. Bouligand, “Sur les singularités ala paroi dans le des ondes liquides,”_ 
= Bulletin des Sciences Mathématiques, France, Vol. 40, 1926, pp. 89-96, 106-112. __ 
Miloh, “Irregularities in Solutions of Nonlinear Wave Diffraction Problem by 7 
™ a Vertical Cylinder,” * Journal of the Waterway, Port, Coastal and Ocean Division, ASCE, — 
‘12. Hunt, J. N., and Baddour, R. E., ‘‘Second-Order Standing Waves Bounded by Circular — 
an Cylinders, ” Journal of the Waterway, Port, Coastal and Ocean Division, ASCE, Vol. 
106, Proc. Paper 15166, Feb., 1980, pp.122-127, 
va 13. Hunt, J. N., and Baddour, R. E., ‘“‘Nonlinear Standing Waves Bounded by Cylinders,” 
Journal of Mechanics and Applied Mathematics, London, England, 4 
33, No. 3, 1980, pp. 357-371. 


Errata. —The following correction should be made to the original paper: 


Page: 290, paragraph I, line Should read “be an imposed pressure 


and z z= = Zr, 0), r > a, the e equation” instead of “be the equation” a 
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PROPULSION OF SPAR Buoys" 


by R. M. Carson 


i _ of experiments using model spar babys have been conducted in 
the regular waves of a wave tank at the Institute of Oceanographic Sciences. 


The buoys often showed a velocity significantly different from that of a small : 


surface marker, or of a dye streak. In many cases the buoy moved faster 
than the surface marker, but sometimes slower, and on occasion it actually 
= . moved toward the wavemaker, against the direction of wave propagation. “The 
internal currents in the tank were not responsible; dye streaks in the water 
showed that the return circulation which wave action in a tank 
was well below the buoy at all times. : | 
_ The explanation of these apparently paradoxical drift velocities seems to lie | 
i the dynamics of the buoy relative to the water surface; these effects are 
over and above the effect detailed in the original paper. ‘Both heave and pitch 
3 The buoy motion in regular wa waves can be characterised by w, and w,, the 
natural frequencies of the buoy in heave and pitch respectively, and by the 
- wave excitation w. In Fig. 3 the expected buoy response is shown diagramatically 


for some alternative combinations of these are 


of course purely schematic. = 


Two main consequencesemerge: 
August, 1980, by Thomas A. McClimans (Proc. Paper 15592). 
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is increased 
of in the buoy velocity | at the crest, and an increase e at the teal. In turn ‘this ae 
will reduce the effect described by the author, 
With the exception of case A, the a vertical component 
of velocity relative to the water; this is especially marked when w > @,- Further, 
the ie buoy may be ata large pitch angle to the vertical when the relative velocity a 
: is a maximum. In this case the relative flow will result in a horizontal component 
of lift force of a By inspection this 
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DISCUSSION” 


not permit a proper of these ideas; 
however, one case is suggestive. A pair of identical buoys were ballasted to 
the same value of draft, but were given different values of pitch inertia. When 
these models were tested side by side in the same regular wave field, it was 

_ found possible to make first one, and then the other proceed faster along the 3 


7 i erser depending on on the wave frequency chosen. This = only be due to their 


_ A further trial was s performed at sea, using spar buoy drifters in a current 


“meter evaluation (4). Floats of markedly different pitch characteristic showed 
very little separation when tracked together over periods of one and of two 
tidal cycles, which Suggests that that wave propulsion was not significant. This negative | 
> result may show that the resonance effects which are all too obvious in the 
regular wave train of a tank become less important in a real seaway. This 
may be due to the fairly sharp resonance (high Q) of the buoy response. —— 


Collar, G., “‘Near- Surface Current Measurement from a "Surface Following Data 
(DB 1)—I,” Engineering, Vol. 5, 1978, pp. 181-196. 


LONGSHORE CURRENTS. Wave C INTERACTION* 


The author is to be complimented for an ‘interesting paper. ‘However, the 

writer would like to clarify a few errors and extend the study to pve 


‘The unsteady longshore momentum equation, in the absence of body Pon 7 


can be written as 93 @ works dole ‘ale — 6 nit 


August, 1980, by Robert A. Dalrymple (Proc. Paper 15592), ts 
Research Engr. _Wave 10. Div. States WES, 


— ud) = (psn) Tox .+—(S,,) 


force is in the directiot™® Popa against the wave 

# 

| 

hout body forces, 

q 


— 
in which u = onshore- velocity. "Using the s definition 
for axes, 1,: 3.6 and 7, and the 


do cos@ | sin@ dx dx 
ak dx 


fort 


-= V(1 + sin? —sin’@ c 


4 au 7 


This equation can be solved analytically to predict the variation of longshore — 
velocity with time. Note that for —" flow sala 0, and if t if the he trignometrical — 


are in ‘Eq. 30, the later reduces | to the author’s Eq. 1. 9 “except the term 
3 sin 6 — 6 sin’ 6, which is an error in Eq. 9, must read 3 sin @ — 4 sin a 
6. In addition if we consider bottom m, to be constant, and neglect 
wave setup or wave set-down, 


alta 


| 
) 
Ez 
» rkC did) 
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ete 
cos 6(1 + sin* 8 


Eq. 35 represents the distribution of aie a currents in the surf zone for 
vt a con constant beach beach | slope. L Likewise carring out th the @ derivatives of 3 of Eq. 23 and 


won, drop sushi 


the 


AIG. 3. —Distribution of Longshore Current (1 ft =0 3048 m) 


Egs. 25-29 into into Eq. 23 23, the unsteady onshore- “offshore momentum 


of 


_ Using the boundary conditions at x = x,, V = V,, in which x, = width 

outer limit of the 

— 


flow a and | since ‘sin 12 = ae sin cos 8, Eq. 37 re reduces 

to the author’s Eq. 8, except in that the term (k?d/4) sin’ @ (dv/dx) — 

be negative, which is an error. Eq. 37 can predict analytically the variation 

In addition, if we consider steady flow, , neglect : wave-setup or wave set- -down, 

a use ~~ 33, ig 37 reduces to a Mater differential equation of the form © 


Xp, V= U,, in which 


QnMx? sin?@ sin” *0 


A of Eqs. 35 40 for k= 0. 78, = des 0.02, 
f = 0.018375, V, = 3.0 fps and U, = 0.3 fps, are shown in oy In addition ‘ 
_ Fig. 3 shows the longshore velocity in to ‘Komar 
~ (10) for the same k, 8, m, f, and V,. 
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WaTER Levet \ VARIATIONS ALONG CALIFORNIA Coast® 


data ia in the or original paper seem to confirm the conclusion reached 
ot by Douglas W. Johnson (13) many years ago, namely that the elevation of 


mean sea level at any location upon the of the connection 


* O’Brien Hall, Univ. of California, Berkeley, Calif. 94720. 
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August, U, by Kaymond A mith and Kobert eltler (Proc. Paper 15014). 
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between the gage and the open « ‘ocean, eee rath ‘the ‘effect of the i i 


hydrography may vary with time due to natural o or man-made changes. 
: _ The paper does not specify the precise location of the tide gages constieent, 
4 but the general character of their exposure to the openocean tide is known. 

In Fig. 5, the gage wk San Diego, located inside the entrance, shows a rise — 
- in sea level with time, whereas the gages at Avila, Santa Monica, and oe 
Angeles, located in or near “open ocean conditions, ‘show: a “nearly | constant 
level. The gage at Crescent City, representing nearly open ocean conditions, — 
_ Showed no change with time whereas the gage at San Francisco, located at 
the throat of the Golden Gate, showed a steady rise withtime. 8” 
* ae The hydraulic regimen, which determines the water surface elevation at the 
throat of an entrance, may be represented with simplifying assumptions by 


During the flood phase, AH is the drop in surface elevation from the ru 
to the throat, sod is the velocity through the throat section, and K is the 
corresponding friction coefficient; on the ebb phase, AH is the drop between &§ 
the bay surface and the throat. The friction coefficients are not neceseasity . 
~ equal for the ebb and flood phases. The surface area of the > bay yy 
flow cross section at the throat Me are variable with h, the | elevation. of the 7 
_ Without examining the hydraulic relationships in more detail, it seems elk 
the difference in mean level between the open ocean and the throat of an_ 
aa is: (1) Not obviously zero; and (2) depends on variables which may 
_ _ Considering the tide-gage record at Fort Point, San Francisco, “the surface 
area of San Francisco Bay has been reduced since 1855, possibly by: as much 
as 30%. The cross-sectional area of this entrance falls near the area-prism curve 
_ of O’Brien (14) and is in equilibrium, the ratio A ,/A . should then be constant, 
but the time required for A, to reach equilibrinm after a reduction in A, is 
not known. The surface width is fixed and a reduction in A. would require 
cana in depth which in turn would increase K. Iti is not difficult to visualize 
_ changes in the surface area and the channels of San Francisco Bay which would 
~ account for the observed increase in the mean level asin Fig.5. = 


_ These notes are an inadequate treatment of a subject which deserves much — ; 


‘theoretical and experimental study. The temporal trend, and the causes, of 


the data available at are » ambiguous. One clear 1 requirement : is an increase in n the 


number of open-ocean tide gages. — 


Johnson, D. W., ‘‘Studies of Mean Sea-Level,” Bulletin of the National Research 
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MIGRATION DUE TO BASELEVEL LOWERING” 


, “Equilibrium Flow Areas of Tidal Inlets on Sandy Coasts,’ * Coastal ig 


Discussion by Saburo Komura, ASCE ne 


: ‘The authors’ extension of the application of wey heat diffusion equation to 
_ a head-cut process is very interesting, but practical application is limited because 
it requires the user to estimate the so-called ‘“‘diffusion coefficient.’’ In — q 
coefficient is unknown and depends on the characteristics of water 
flow, sediment transport and s so on. writer believes. that in the ‘original 
paper there is considerable room for improvement in the theoretical Gevelogmemt — 
, _ of very complex head cut process. In this discussion, the writer will concentrate a 
on the determination method of diffusion coefficient for the head-cut migration © 


Bes Substituting Eq. 17 into > Eq. 15, the following equation of motion of a headcut 


ex| -(; =) | — ad 


ex 19) 


in which w = « dx/dt (propagation v velocity < of headcut); j = -Y,/h which is 
determined by the outlet conditions; and b = the initial bed slope. The writer 
believes the second term of Eq. 18 should read bk/h instead of bk. From — 
19, the following dimensionless can obtained: tind 
ont fal 


\ 


Taking 1g Vkt/h as 
be expressed as shown in Fig. 19. The equation of envelope appeared in Fig. A 
19 was obtained from the following two equations: qint dee 


}=0 


= “August, 1980, by Ze’ ev B B. Begin, David F. Meyer, and Stanley A. Schumm (Proc. 


4 
of | 
| 
&g | 
=] 
q 
4 
4 
| 


Substituting 4 ‘nin into 20a 


© 
which e= ‘the of natural Also, substituting Eq. 


(7 


=" 


‘anal ty 


- 


4, 


= which Y = V kt/h. Differentiating Eq. 20a with respect to Y, the pard a 
— |—] +5). (25a) 
| 
\ 
b=0.015 
— ° 


~ Observation results on head-cut migrations can be oieiaenl in athe — 


"Propagation velocity is obtained by differentiating Eq. 26 with respect 


a 


By 26 and 27, Eq. 19 becomes 


The k should be time- independent, for small value of the 
_k-value obtained from this equation will be time- dependent. On the other hand, “Ts ie 


for B = | and a large value of ¢ 
Accordingly, the following equation is obtained: 
When aie x axis is re along the i ‘initial bed surface, or b=0 (very small ad 


3 = 1/2. 


centimeters | centimeters | centimeters | centimeters a centimeters. 
per minute | per minute | per minute per minute per minute 


23.4(0.69) 
7(0. 
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ISCUSSION 


‘The ‘authors’ data were analyzed by using Eqs. 31 and 33, and data for b ae 
0.015 (Run 4, 5, and 12, and B = 1) were plotted in Fig. 19 ie ee 
a Finally, the writer would like to mention the terminal time and distance oc 
headcut. Assuming that wh/k at the terminal time has a value to 


mek 


a value of j Pp pone be 10% | OF So. Also, from Eq. 25a, the terminal distance 


i 
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RAT! AVES AN WaAvE Turorres* 
Discussion by D. H. 


se water by a pneumatic wave generator and allowed to propagate over a 
region of variable depth onto a region of constant depth where the measurements 
ae were made. Since nothing is stated about their generation, ‘it seems likely that — 
simple sinusoidal air pressure was applied. ‘It is clear, e.g., , from Figs. 
ii 5, 8(a), 8(b), 8(c) that the waves do not correspond to uniformly propagating — 
waves and thus change in profile along the region of constant depth. A more 


a4 valuable ccmparison with the theories for regular periodic waves would na? 
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8 q 
Given values of j = 
Ps and error. The initial value of k for a trial-and-error method can be estimated oe | 
5 @6€=—Ss—téSTThhiiss equation was obtained from Eqs. 22 and 26, i.e., x = V 2kt = at®, since — | 
«B= 1/2. Table 1 shows values of a, B and k. It would be interesting to note 
. 


been obtained if some to make them been ‘betta “ea 
see Hansen and Svendsen (10). Otherwise, as the author states in his ‘inal 
: 1 sentence, the data refers only | to the particular configuration of that tank. ; 
: _ The comparisons with theory clearly show that the author has failed to realize 
that in order to specify a wave solution completely the mean flow and mean 
a depth must be determined as well as a suitable set of wave parameters. The 
difficulty this neglect has lead to is clearly evident in the comparisons of the 
profiles of uw, and in Figs. 8(b), and 8(c) where it is clear that 
_the * ‘theoretical’ lines all correspond to different values of these mean quantities. 
The importance of clearly defining the mean flow, or equivalently the reference 
; frame, for which theoretical work is developed has been recognized since Stokes 
3 drew attention to it. More recently Jonsson (11) has drawn attention to” 
: its importance. The set-up and set-down due to water waves propagating from 
- depth to another is also a a well well known phenomenon, e.g., see. Bowen, et 
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